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I.  INTRODUCTION 

A  new  class  of  integrated  circuits,  called  charge  domain  devices  (C'I)D)  has  been  developed  with  the 
goal  of  performing  signal  processing  functions  with  accuracy  anil  speed  performance  exceeding  aherna 
live  technologies.  The  starting  point  for  this  development  is  conventional  charge  transfer  ilevue 
(CIO)  technology.'1'  Devices  of  this  type  such  ax  charge  coupled  device  K CD)  transversal  filters' 
have  been  demonstrated  in  many  cases,  particularly  at  high  frequencies,  to  he  mote  ellicieni  m  pci 
forming  particular  sampled  data  processing  functions  than  such  alternative  techniques  as  digital  liltcrx  m 
switched  capacitor  devices.  However,  as  the  speed  and  accuracy  performance  requirements  increase, 
conventional  CTO  s  also  encounter  a  mber  of  limitations,  and  it  is  the  purpose  of  the  development 
effort  covered  in  this  report  to  provide  a  technology  which  permits  higher  speed  operation  and  wider 
dynamic  range  by  avoiding  these  limitations.  In  order  to  understand  the  origin  of  some  of  these  limita¬ 
tions.  the  operation  of  conventional  CTOs  are  briefly  reviewed. 

Conventional  CTO's  derive  their  output  signal  by  sensing  charge  packets  with  overlying  MOS  elci- 
trodes.  Multiplication  of  these  charges  by  tap  weights  is  implemented  by  splitting  the  overlying  elec¬ 
trodes  in  proportion  to  the  desired  impulse  response  coellicienl,  and  summation  is  implemented  In  con¬ 
necting  the  overlying  electrodes.  This  implementation  is  particularly  ellicieni.  since  all  of  the 
mathematical  operations  are  accomplished  automatically  by  simple  physical  laws  ralhci  than  by  manipu¬ 
lating  binary  bits  in  complex  logic  circuits,  but  it  causes  at  least  two  of  the  problems,  first,  as  the 
specifications  on  the  system  increase,  and  the  number  of  filter  coefficients  needed  to  accomplish  the 
desired  transfer  function  increases,  the  total  capacitance  of  the  output  electrode  increases,  making  high 
speed  operation  more  difficult;  and  second,  non-linearities  in  the  relationship  between  the  charge  in  the 
packets  and  the  voltage  induced  on  the  overlying  electrodes  usually  compromises  the  accuracy  of  the 
transfer  function  if  buried  channel  technology  is  used.  These  constraints  have  in  the  past  limited  the 
frequency  handling  capabilities  of  conventional  CTO  filters  to  a  few  megahert/. 

Charge  domain  integrated  circuits  have  been  developed  to  overcome  these  limitations  so  as  to  in¬ 
crease  the  frequency  range  that  can  be  handled  by  monolithic  signal  processing  chips.14  '1  In  charge 
domain  devices,  all  signal  processing  is  performed  by  manipulating  the  charge  packets  themselves,  ratti¬ 
er  than  using  the  image  charge  on  overlying  electrodes.  The  charge  packets  representing  the  input  sig¬ 
nal  may  be  split,  routed,  delayed  and  combined  to  form  new  charge  packets  that  represent  the  output 
signal.  But  since  the  splitting  and  routing  depend  only  on  the  plan  view  geometry  of  the  devices  that  ac¬ 
complish  it,  and  not  on  the  details  of  the  capacitance-voltage  characteristic,  buried  channel  technology 
can  be  used  without  degrading  the  accuracy  of  the  transfer  function,  furthermore,  since  portions  of 
some  charge  packets  can  be  routed  backwards  in  the  signal  flow  sense  and  re-introduced  into  the  for¬ 
ward  path,  CDDs  yield  the  new  possibility  of  filters  with  infinite  impulse  response;  i.e„  tillers  that  im¬ 
plement  poles  as  well  as  zeroes  in  their  transfer  function,  finally,  since  the  output  of  these  devices  is  a 
stream  of  charge  packets,  low  capacitance  diode  sensing  of  the  charge  is  used  to  generate  the  output  sig¬ 
nal.  Thus,  the  output  capacitance  does  not  increase  as  the  tiller  architecture  becomes  more  complex, 
and  dev  ice  speed  is  limited  only  by  the  speed  of  charge  transfer  which,  for  buried  channel  technology, 
may  be  as  large  as  hundreds  of  megahertz. 

The  purpose  of  this  contractual  effort  was  to  investigate  the  unique  capabilities  of  charge  domain  de¬ 
vices  and  to  develop  them  for  high  speed  sampled  data  analog  signal  processing.  Specifically  the  three 
goals  of  the  contract  tire,  as  slated  in  the  original  proposal:"’1 

•  To  investigate  the  fundamental  physical  operation  and  limitations  of  the  charge  domain  concept  and 
the  relationship  to  analog  tiller  performance  characteristics.  This  will  cover  the  relationship  between 
such  device  characteristics  as  lateral  transfer  speed,  charge  packet  splitting  accuracy,  and  such  per¬ 
formance  characteristics  as  filter  accuracy,  linearity,  and  bandwidth. 

•  To  design  and  implement  a  charge  domain  filler  as  a  lest  vehicle  that  will  demonstrate  performance 
capabilities  loi  a  selected  application  I  he  filler's  design  will  include  the  lull  laupe  ol  liade  oils  be 
tween  poles  and  /cuts,  a  unique  fe.iluie  ol  the  i  h.upe  domain  concept 

•  lo  investigate  the  feasibility  of  iniioducing  programmability  into  charge  domain  Idlers  Basu  pin- 
grammable  building  blocks  will  be  designed,  fabricated,  anti  evaluated  as  part  of  this  program.  It  is 
envisioned  that  the  results  of  this  work  will  lead  lo  a  future  program  to  pioducc  a  digitalli  piopi.un 


mahlc  Idler  module  contiollcd  h>  readily  available  micrnprocexxois  In  iliis  fashion.  Inline  signal 

processing  building  blocks  can  ho  implemented  will)  software  configurability  anil  thus  quickly  ami 

conveniently  integrated  into  a  variety  of  applications. 

This  investigation  has  covered  all  the  areas  outlined  in  the  proposal  and  has  gone  inio  more  depth  in 
several  of  the  areas.  The  goals  of  the  contract  have  been  addressed,  and  the  successful  demonstration 
of  design  methodologies,  test  vehicle  filters  and  programmability  concepts  has  proven  that  the  charge 
domain  concept  is  indeed  a  powerful  new  technique  for  flexible  realizing  signal  processing  functions, 
and  ultimately  entire  systems  on  a  single  integrated  circuit  chip. 

Section  2  discusses  some  of  the  basic  charge  domain  operations  for  delay,  multiplication  and  addi¬ 
tion.  Section  2  describes  the  representation  of  signals  in  (he  charge  domain  which  permits  four  quad¬ 
rant  arithmetic  functions  to  be  reali/ed.  Also  in  this  section  is  a  description  of  the  mathematics  under¬ 
lying  a  filter  design.  This  leads  into  Section  4  which  details  some  of  the  device  structures  which  are 
needed  to  implement  the  mathematical  functions  and  their  capabilities  and  limitations.  Several  charge 
domain  structures  which  overcome  these  limitations  are  described  in  Section  5. 

Section  6  discusses  performance  characteristics  such  as  linearity,  dynamic  range  and  bandwidth.  The 
operation  of  charge  splitting,  which  determines  coefficient  accuracy  is  also  discussed.  The  physics  ol 
this  operation  has  been  investigated  through  testing  of  structures.  These  results  indicate  experimentally 
how  the  splitting  structure  may  be  optimized  for  accuracy. 

Section  7  describes  computerized  test  stations  which  have  been  developed  for  data  acquisition  to 
characterize  the  accuracy  and  noise  levels  of  charge  domain  devices.  Experimental  results,  which  are 
much  more  precise  than  those  obtained  by  previous  methods,  are  also  presented. 

Section  X  discusses  the  computer  aided  filter  design  tools  which  have  been  developed  as  the  work 
proceeded.  These  tools  tire  used  in  the  interpretation  of  the  test  results  and  also  in  determining  the 
feasibility  of  certain  filter  designs  prior  to  device  fabrication. 

Section  4  details  the  actual  design  and  operation  of  two  charge  domain  filters  which  have  been  de¬ 
veloped  for  this  contract,  a  simple  three  pole  low  pass  filter,  and  an  eight  pole  narrowband  bandpass 
filter,  limit  of  these  fillers  have  been  successfully  demonstrated  and  experimental  measurements  of 
their  performance  are  included. 

\s  mentioned  above,  one  of  the  major  goals  of  this  contract  is  to  provide  the  capability  of  electrical¬ 
ly  programming  the  coefficients  that  characterize  the  transfer  function  of  a  C'DD  filler.  Section  10  de¬ 
scribes  several  approaches  which  can  be  taken  to  realize  digitally  programmable  devices.  Two  of  these, 
namely  a  sequential  method  which  minimizes  chip  area  and  a  pipelined  method  which  optimizes  speed, 
have  been  successfully  implemented,  and  experimental  results  are  presented. 

Section  I  I  looks  to  the  future  and  describes  other  features  of  the  charge  domain  which  have  not  yet 
been  utilized.  I  Ins  section  describes  some  applications  which  could  be  efficiently  addressed  with  charge 
domain  technology  and  would  eventually  result  in  the  desired  system  on  a  chip.  All  of  the  results  to 
date,  and  possible  future  extensions  are  summarized  in  Section  12. 
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2.  BASIC  CHARGE  DOMAIN  OPERATIONS 


Discrete  time  signal  processing  requires  only  a  relatively  small  number  of  functions,  namely,  dela\ 
(memory),  and  four  quadrant  addition  and  multiplication.  In  order  to  opiimi/c  performance,  these 
functions  must  be  implemented  over  a  wide  dynamic  range,  yet  consume  little  power  and  silicon  area 
This  section  describes  how  the  charge  domain  concept  addresses  these  requirements 

2.1  INTACT  CHARGE  PACKET  TRANSFER 

Charge  packet  transfer  is  the  basic  operation  used  in  both  the  CCD  delay  line  ami  the  CCD  spin 
electrode  transversal  filter.  It  can  be  accomplished  in  a  number  of  different  wavs.' ' M  flic  technique 
used  for  this  contract  is  the  pseudo  one  phase  approach  shown  in  top-view  and  cross  section  in  I  ie- 
ure  2.1  This  technique  requires  4  electrodes  per  stage.  A  50%  duly  cycle  clock  is  applied  to  the  two 
electrodes  in  the  center  of  figure  2.1.  These  electrodes  are  called,  respective!) ,  the  clocked  iranslci 
electrode  and  the  clocked  storage  electrode.  The  other  two  electrodes  are  held  ai  DC.  and  are  called 
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f  igure  2.1  Pseudo  one  phase  CCD  operation 
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I"' I  1 1 1 1 1 1  "II  Hi',  tii.ni'v.  |i.nkd  iiaiistcis  in  die  light  underneath  die  IK  sioiagc  I'la'inilc  I ‘ms 
i lie  10  is  one  lull  stage  ol  iransler  per  one  lull  clock  cycle  (one  delay  period)  As  a  shorthand  notation 
lor  other  figures.  one  lull  stage  including  four  electrodes  is  shown  schematically  as  a  single  rectangle 
representing  the  plan  view  of  the  reservoir  holding  the  charge  during  one  of  the  clock  phases 

2.2  PAt  KKI  SPUTUM, 

\s  si, own  in  figure  2.2.  this  operation  is  identical  lo  the  basic  operation  described  above  evvepi  ilia, 
i lie  receiving  reservoir  is  comprised  ol  two  or  more  sections  whose  total  capacitance  is  roughly  equal  m 
dial  ol  the  source  reservoir,  but  whose  relative  si/es  correspond  to  predetermined  i.iiios  <>mc  iiv 
charge  lias  been  split  into  portions,  each  individual  portion  can  be  mampulaied  unlepeinlenilv  I  no 
unilaieral  splitting  action  can  be  conveniently  obtained  by  starting  the  leading  edges  ol  i lie  bainet  ic 
gioiis  that  separate  l he  receiving  reservoirs  in  the  transfer  gale  region  of  die  slnlt  registei  struvtuic  as 
shown  in  the  figure  If  die  potential  in  the  transfer  gale  region  is  always  higher  Ilian  die  icsciving 
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f  igure  2.2.  Barrier  charge  packet  splitting. 
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reservoirs,  the  charge  will  nol  he  able  lo  “hack  up"  around  the  harriers  between  the  reservoirs. 
Mlhough  the  leading  edges  of  the  harriers  are  shown  coming  to  a  point,  tins  feature  is  not  essential  li 
is  helplul.  however,  in  minimizing  the  effect  of  small,  uneontrolled  electric  fields  on  the  division  of  tire 
loial  charge  among  the  receiving  reservoirs. 


I  he  packet  splitting  operation  described  above  is  used  lo  implement  multiplication,  but  some  prob¬ 
lems  must  he  laced  First  of  all,  it  is  not  be  possible  to  implement  coefficients  greater  than  unity  with 
th's  approach,  and  second,  the  implementation  of  multiplication  by  a  negative  number  is  not  obvious, 
both  ol  these  points  are  addressed  in  detail  later. 


2. A  I’M  kill  COMBINING 


\s  win  he  seen  in  Figure  2  A.  this  operation  is  also  essentially  identical  to  the  basic  transfer  opeta- 
tion  ihskUsseil  above  I  his  lime  several  component  packets  are  simultaneously  transferred  into  a  single 
reset v on.  and  n  is  obvious  that  this  operation  implements  the  summation  of  the  individual  component 
i  bulge  packets 


2.4  KQl  Il  IBK  Vf ION 


A  technique  called  "equilibration"  is  used  to  implement  both  addition  and  multiplication  in  a  single 
step  In  this  operation,  which  is  shown  in  Figure  2.4.  two  charge  storage  regions  arc  briefly  connected 
together  during  each  clock  period  by  turning  on  the  equilibration  gale  that  separates  them.  During  the 
time  this  gate  is  open,  charge  Hows  between  the  reservoirs  so  as  lo  equalize  (heir  potentials.  After 
equilibration,  the  gate  is  turned  off,  leaving  amounts  of  charge  in  each  reservoir  in  proportion  to  their 
respective  capacitance  (assuming  equal  applied  voltages  and  threshold  voltages)  irrespective  of  the  ini¬ 
tial  distribution  of  charge.  Since  charge  is  required  to  move  in  both  the  A  and  >  directions  in  this 
structure,  there  is  a  potential  speed  problem.  The  relatively  long  time  required  for  an  initially  unequal 
charge  distribution  to  become  uniformly  distributed  can  be  greatly  reduced  by  introducing  a  lateral 
diffusion  channel  to  lower  the  lateral  resistance  and  thereby  reduce  the  time  constant.  This  diffusion 
may  be  overlaved  by  a  metal  to  further  enhance  lateral  transfer  speed.  The  penalty  inherent  with  this 
technique,  which  results  in  a  structure  similar  to  a  stage  in  a  bucket  brigade  delay  line,  is  higher 
transfer  loss  However,  since  only  a  small  number  of  these  stages  are  required  in  a  device,  this  is  usu¬ 
ally  not  a  limitation 
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l  igure  2.3.  Charge  package  summation 
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3.  SIGNAL  PROCESSING 

3.1  SIGNAL  REPRESENTATION 


i or  general  li Iter  design,  it  is  necessary  to  represent  positive  ;md  negative  impulse  response 
coefficients  as  well  as  positive  and  negative  signal  values,  and  it  is  necessary  to  implement  mullipliea 
lion  correctly  in  all  lour  quadrants.  Although  representation  of  positive  and  negative  signal  values  can 
he  accomplished  in  the  charge  domain  by  the  use  of  a  bias  charge,  it  is  mote  convenient  to  use  two  sep¬ 
arate  charge  packets  for  each  signal  sample.  In  this  case,  the  actual  signal  level  can  be  represented  In 
the  dillerence  between  these  two  charge  packets  This  approach  has  the  advantage  that  the  zero  level  is 
not  required  to  be  fixed,  but  can  be  at  any  level;  all  that  is  required  to  represent  a  zero  signal  level  is 
that  the  two  packets  be  equal  to  each  other.  Thus,  all  of  the  signal  processing  elements  which  aie  da 
scribed  have  two  input  charge  packets  for  each  signal  sample,  and  they  will  produce  two  output  eliuige 
packets  I  he  output  signal  sample  is  represented  by  the  dillerence  between  the  charges  in  these  uuipui 
packets  without  regard  to  their  absolute  values.  Thus,  these  output  packets  are  in  the  exact  same  lor- 
mat  as  the  input  packets,  and  can  be  fed  into  another  signal  processing  element  directly  to  implement 
the  cascade  of  the  processes.  A  differential  amplifier  is  required  to  produce  the  final  output  signal,  but 
no  amplifiers  are  needed  at  any  intermediate  point  in  the  structure.  This  representation  provides  the 
additional  benefit  of  rejecting  common  mode  interference  at  the  output  siu.ii  as  feedthrough  ol  clock 
waveforms,  dark  current,  and  other  undesired  signals. 

Packet  splitting  is  used  to  implement  multiplication.  Positive  co-efficients  (less  than  unity)  are  im¬ 
plemented  by  delivering  the  desired  fraction  of  each  packet  to  the  corresponding  output  packet,  and 
negative  coefficients  (less  than  unity)  are  implemented  by  delivering  the  selected  portions  to  the  oppo¬ 
site  polarity  output  packet.  It  is  easy  to  show  that  this  scheme  does  indeed  implement  four  quadrant 
multiplication  correctly  (Note  that  a  negative  signal  is  represented  by  a  larger  negative  charge  packet 
and  that  this  is  delivered  to  the  positive  output  when  the  coefficient  is  negative.  Tint-.,  a  negative  signal 
multiplied  by  a  negative  coefficient  results  in  a  positive  output  ) 

3.2  IMPLEMENTATION  OF  MATHEMATICS 

A  general  form  of  the  difference  equation  that  is  to  be  implemented  by  the  filter  is. 
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and  the  corresponding  transfer  function  for  the  /.-transform  is 
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The  poles  of  the  transfer  function  are  the  roots  of  the  polynomial  in  the  denominator,  while  the  /eros 
are  the  roots  of  the  numerator  in  Equation  (3.2). 

To  keep  the  physical  implementation  simple,  we  restrict  the  coefficients  to  real  numbers.  It  is  often 
useful  to  consider  the  denominator  as  a  product  of  first  and  second  order  factors  having  real 
coefficients.  Since  the  original  rational  transfer  function  can  be  implemented  as  a  cascade  of  individual 
factors  or  as  a  sum  of  partial  fractions,  it  is  only  necessary  to  discuss  the  implementation  of  one  first  or¬ 
der  and  one  second  order  factor.  Although  the  numerator  can  also  be  factored  in  a  similar  manner,  it  is 
usually  practical  to  implement  several  zeros  at  once  by  using  transversal  tillers  of  modest  length.  Thus, 
it  will  not  be  necessary  to  deal  with  the  elementary  factors  in  the  numerator.  Also,  note  that  numerator 
coefficients  can  always  be  scaled  such  that  their  sum  is  unity,  and  that  the  only  effect  is  a  change  in  the 
overall  gain.  Unfortunately,  the  same  cannot  be  said  for  the  denominator.  The  unit  term  in  the 
denominator  of  equation  (3.2)  came  from  the  left  hand  side  of  equation  (3.1),  and  this  term  sets  the 
scale  lor  the  ,-f’s.  Furthermore,  most  of  the  desirable  pole  locations  involve  second  order  factors  in 
which  the  coefficient  of  the  first  order  term  has  a  magnitude  greater  than  unity.  Since  it  is  impossible 
to  implement  such  coefficients  in  a  purely  charge  domain  design,  it  has  long  been  thought  that  these 
poles  cannot  be  implemented  without  resorting  to  devices  that  provide  active  gain. 


4.  GENERAL  APPLICATION  REQUIREMENTS 

The  charge  domain  basic  building  block  (unctions  can  be  connected  lust  as  the  sections  of  a  block 
diagram  of  a  sampled  data  filter  are  connected  in  order  to  provide  for  the  poles  and  /ernes  of  the  fillet's 
transfer  function.  There  tire  limitations  on  the  placement  of  the  poles  and  /ernes  due  to  the  constraint 
of  having  no  gain  elements  in  the  system  and  of  performing  only  fractional  multiplications  I  he  real 
axis  pole  corresponds  to  an  exponential  decay  m  the  time  domain,  while  the  complex  pole  pan 
corresponds  to  a  decaying  sine  wave. 


I  he  desired  transfer  function  for  a  real  axis  pole  is 
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and  so  the  pole  is  placed  at  Z  =  ii|.  The  Z  1  term  in  the  numerator  of  liquation  4.1  only  provides  un¬ 
necessary  delay  of  the  output.  Removal  of  this  delay  by  placing  a  zero  at  the  origin  <Z  =  0)  xx ill  pro¬ 
vide  a  smaller  CDF  structure  and  reduce  the  throughput  delay.  The  resulting  dillerence  equation  is 


HZ)  =  <»,Z  1  HZ)  +  A,.V(Z)  . 


(4.2 ) 


The  block  diagram  for  this  function  is  shown  in  Figure  4.1.  By  varying  «|  from  zero  to  one  or  zero  to 
minus  one,  the  pole  can  be  placed  anywhere  on  the  real  axis  within  the  unit  circle  (-1  <  Z  <  I).  At 
the  split  the  sum  of  the  magnitudes  of  the  coefficients  must  be  less  than  or  equal  to  one  so  the  max¬ 
imum  gain  is  when  A  =  1  -  l«|l.  This  gain  constraint  can  be  compensated  for  by  proper  sizing  of  the 
adjacent  circuitry.  Thus  implementation  of  a  pole  on  the  real  axis  is  realizable  within  charge  domain 
constraints. 


The  transfer  function  of  the  complex  pole  pair  is 


H(z) 


ki _ 

z~  —  2a2z  +  (a2  +  ) 


_ k{r2 _ 

1—2 a2Z'  +  («?2  +  wj ) Z~ 2 


(4.3) 


thus  placing  poles  at  Z  =  a2  ±  jw2-  The  Z ~2  term  in  the  numerator  can  be  removed  with  a  second  or¬ 
der  zero  at  the  origin.  The  resulting  difference  equation  is 


HZ)  =  A  >  V  ( Z )  t  2o,Z  'HZ)  -  (or  +  h  /  >Z  ’HZ) 


(4  4) 


Its  block  diagram  is  given  in  Figure  4.2.  In  this  case  the  sum  of  the  required  coefficient  magnitudes  at 
the  split  may  become  greater  than  one  even  for  A ,  =  0  for  certain  pole  locations.  The  realizable  pole 
locations  for  direct  charge  domain  implementation  of  a  second  order  section  are  inside  of  the  dashed 
lies  of  Figure  4.3.  This  imposes  a  severe  limit  on  design  flexibility  since  many  of  the  desired  poles  are 
outside  of  this  region. 

fhe  solution  i s  to  lind  functions  which  can  be  implemented  in  t he  charge  domain  which  provide 
poles  m  their  desired  locations  plus  some  other  non-desired  poles.  Ihe  spurious  poles  may  he  cancelled 
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I  igtire  4  1  Recursive  accumulaloi  block  diagram 


out  with  /crocs  which  can  he  placed  anywhere  in  the  /-plane  by  using  a  charge  domain  transversal  filler . 
One  such  function  is 
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The  associated  difference  equation  (assuming  the  addition  of  an  Nth  order  /ero  at  the  origin)  is 


HZ)  -  k;X(Z)  -  AZ  '  HZ) 
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where  A2  =  1~A.  The  block  diagram  of  a  recursive  resonator  which  implements  this  function  in  shown 
in  f  igure  4.4.  The  roots  of  this  equation  are 

Z  =  (  a)"  '  ii  -  1 ,2,  •  •  •  V  (4  ’> 


which  yield  N  poles  equally  spaced  around  the  unit  circle.  In  polar  coordinates  the  poles  are  located  .< 

Z  =  o'  exp  \H2n-  \)tt/2N]  US) 

for  example,  with  ,V  =  5,  A  =  0.2  and  A_>  =  0.8,  the  poles  are  located  as  shown  in  figure  4.5.  If  one  is 
allowed  to  freely  choose  N  and  A,  poles  can  be  placed  anywhere  within  the  unit  circle.  Note  that  it  is 
impossible  to  build  an  unstable  filler  in  the  charge  domain  because  the  magnitude  of  the  coefficient  >>  is 
limited  to  less  than  one  by  lack  of  gain  elements  in  the  circuit  and  the  conservation  of  charge  at  the 
splitting  junction.  Normally,  one  is  interested  only  in  a  small  number  of  the  pole  pairs  produced  by  the 
structure  shown  above,  and  the  unwanted  poles  must  be  canceled  by  a  transversal  tiller.  When  the 
zeroes  implemented  by  the  transversal  filter  are  placed  directly  on  top  of  the  unwanted  poles,  the 
desired  transfer  function  is  obtained. 


figure  4.4.  Recursive  resonator  block  diagram. 


5.  ('HAR<;K  DOMAIN  RKCURSIVK  STRUC TURKS 


A  structure  which  implements  the  re;il  pole  function  is  shown  schematically  in  figure  Vl  I  lie  im¬ 
pulse  or  unit  sample  response  should  be  a  decaying  exponential  I  he  charge  packet  representing  .1  unit 
sample.  \.  is  clocked  into  Section  A.  I  he  equilibration  gale  opens  and  closes  leaving  H  with  ,,  \  .uni 
I  w  ith  A ,  A  .  w  here 

O  |  -  (  /if  I  (  |  t  (  /( )  ( s  |  1 

and 

A  1  —  l—o  |  <  1/  <  (  ,  f  (  if).  ('At 

The  portion  under  I  is  clocked  out  and  no  new  charge  enters  (for  a  unit  pulse  input).  I  lie  equilibra¬ 
tion  gale  opens  and  closes  again  storing  o  f  A  in  H  and  outputting  <1 ,  A  ,  \  This  pmcevlure  continues  lm\ 
mg  the  required  exponential  decay  of 

1 1  n  1  A  a-  A  1  s  ; 1 

where  1  (/;)  is  the  output  sample  and  n  is  the  sample  number. 

A  charge  domain  structure  called  a  circular  shift  register  was  developed  prior  to  this  contract  to  real¬ 
ize  the  transfer  function  of  liquation  4.6. <41  This  structure  is  shown  .schematically  in  f  igure  v2  Since 
it  is  to  be  a  member  of  a  family  of  devices  in  which  signal  samples  are  uniformly  represented  as  the 
dillerence  between  two  charge  packets,  it  is  shown  with  two  input  structures.  Ibis  is  a  key  feature 
which  is  missing  from  similar  approaches  developed  by  others."”  In  the  figure,  the  positive  weight 
charge  packets  enter  at  the  lop  from  the  left,  and  the  negative  weight  charge  packets  enter  from  the  left 
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Figure  5.2.  Circular  shill  register  schcmalic  showing  impulse  response. 

al  the  hotiom.  Similarly,  the  positive  weight  output  packets  leave  at  the  top  toward  the  right,  and  the 
negative  weight  packets  leave  toward  the  right  at  the  bottom. 

I  <>ur  ipi.idrant  operation  ti  e.,  both  positive  anil  negative  inputs  and  coefficients)  is  implemented  by 
using  this  dillerenlial  structure  and  by  representing  a  positive  signal  sample  by  a  charge  packet  ol 
greater  than  the  50%  ol  saturation  level  on  the  positive  side  and  a  corresponding  charge  packet  on  the 
negative  side  having  less  than  50%  amplitude,  while  a  negative  signal  is  represented  by  a  positive 
weight  charge  packet  having  less  Ilian  50%  ol  the  lull  potential  well  and  a  negative  weight  packet  of 
mine  than  amplitude  I  he  linal  output  will  be  formed  by  taking  the  dillerence  between  the  two 
output  packets.  Hi)  )  t  >.  which  also  provides  enhanced  common  mode  rejection.  The  operation 
ol  the  v iiciilar  shift  registei  may  best  be  understood  by  looking  at  Us  impulse  response,  i.e..  tracking  the 
output  vcisiis  tunc  tin  clock  periods)  for  a  single  full-scale  input  pulse  at  the  V<  +  )  input  at  lime  sam¬ 
ple  a  <1  and  will)  all  other  inputs  zero.  I  he  impulse  charge  packet  proceeds  to  point  I  where  it  ts  di¬ 
vided  In  a  barrier  splitter  with  a  portion  I  -  A  proceeding  into  the  positive  output  channel  and  Ihe  A 
portion  following  the  feedback  channel  I  lie  first  output  pulse  appears  with  tin  amplitude  of  I  -  A 
I  ive  cycles  later  the  A  portion  litis  proceeded  to  point  H  where  it  is  again  split  into  A  and  I  -  A  portions. 
The  output  portion  appears  as  -  Ad  A)  with  \  '  traveling  back  up  towards  point  A.  Five  cycles  later 
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Figure  5.3.  Circular  shift  register  impulse  and  frequency  responses. 

the  split  occurs  at  A  once  again  with  a  resulting  output  of  \2(1-X),  and  so  on.  This  impulse  response 
(which  is  the  inverse  Z  transform  of  Equation  4.6)  is  shown  in  Figure  5.3.  Note  that  any  periodic  sig¬ 
nal  whose  period  matches  the  time  needed  to  complete  one  full  circle,  and  which  changes  sign  during 


each  half  period  adds  onto  itself  as  it  circles  around  the  shill  register.  Thus,  a  sine  wave  with  a  period 


equal  to  the  number  of  stages  around  the  circular  shift  register  builds  in  amplitude,  and  the  structure 
also  resonates  at  all  odd  multiples  of  this  fundamental  resonance.  These  additional  harmonic  reso¬ 
nances  correspond  to  the  additional  “unwanted”  poles.  All  other  frequencies  do  not  resonate  and  are 
attenuated.  The  resulting  frequency  response  is  also  shown  in  Figure  5.3. 

A  charge  domain  transversal  filter  can  be  cascaded  with  the  circular  shift  register  to  provide  zeroes 
to  cancel  out  all  the  higher  harmonics  leaving  the  fundamental,  or  cancelling  out  the  fundamental  and 
leaving  one  of  (he  harmonics.  The  form  of  a  charge  domain  transversal  filler  for  leaving  only  the  fun¬ 
damental  is  shown  in  Figure  5.4.  Its  impulse  response  is  a  sampled  data  vcision  of  one  half  cycle  of  a 
decaying  sine  wave.  When  this  impulse  response  is  convoked  with  that  ol  the  circular  shift  register, 
the  oiled  is  to  till  in  the  empty  samples  of  Figure  5.3  with  samples  ol  the  desired  waveform  I  he  im¬ 
pulse  response  of  this  cascaded  system  is  the  sampled  data  decaying  sine  wave  shown  in  I  igtire  5.4 
This  truly  is  an  infinite  impulse  response  design.  The  resulting  frequency  response  is  shown  in  Fig¬ 
ure  5  5.  Note  that  now  only  the  resonance  at  the  fundamental  frequency  is  apparent,  till  other  frequen¬ 
cies  are  attenuated.  I  his  structure  is  capable  of  implementing  poles  essentially  anywhere  within  the 
unit  circle  ol  the  /-plane,  with  the  angles  expressible  as 

«  -  Imi/m  (s.4l 


s* 


where  n  and  in  are  integers 

While  this  technique  places  poles  anywhere  within  the  unit  circle,  in  some  cases  this  approach,  while 
feasible,  is  not  practical.  Suppose  for  instance  a  complex  pole  pair  is  to  be  placed  at  angle  of  ft  =  ±  89" 
in  the  /-plane.  According  to  Equation  4.8  this  requires  a  circular  shift  register  with  2 N  —  180  stages. 
The  transversal  tiller  is  required  to  cancel  out  88  unwanted  poles.  In  addition  to  being  impractical  in 
terms  of  size  and  power,  the  pole  section,  even  a t  its  non-resonance  frequencies  provides  little  attenua¬ 
tion,  and  thus  the  entire  filter  is  only  marginally  better  than  the  transversal  filter  alone. 


A  modification  of  llic  circular  shift  register  concept  has  been  discovered  which  provides  a  much 
more  elficienl  way  to  place  poles,  with  only  a  slight  loss  ol  llevibihlv  A  stand. nd  ciuulai  slid!  ugistet 
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I  ipurc  5  4  Charge  domain  transversal  filter  schematic  and  overall  impulse  response. 

vMilt  2  A  12  Slavics  in  the  feedback  loop  has  ,V  =  6  poles  at  angles  of  n  -  (2//  -  I)  30'',  n  I.  2,  ...6. 
A  standard  circular  shill  regisier  with  2.V  -  14  stages  places  N  —  7  poles  at  angles  of  t<  =  (2m  I1 
2s  7'  ||  ihc  circular  shill  regisier  is  modified  lo  include  multiple  feedback  loops,  the  poles  mat  be 

placed  more  elleciively.  for  instance  some  portion  "A"  may  be  fed  back  with  12  stages  and  some  por¬ 
tion  "  R  led  back  with  14  stages,  where  I  t  R  <  I.  By  varying  the  4  to  R  ratio  the  poles  can  be 
placed  anywhere  on  the  locus  ol  points  between  the  2N  =  12  and  IN  -  14  locations  as  shown  in  fig¬ 
ure  s  (.  I  bus  instead  ol  requiring  180  stages  of  feedback  and  88  cancellation  zeroes,  this  technique 
only  requires  2\  14  singes  will)  5  cancellation  /crocs  to  place  a  pole  at  o  =  89".  The  only  drawback 

is  dial  these  poles  can  no  longer  be  placed  arbitrarily  close  to  the  unit  circle  even  for  A  t  R  ^  I. 

I  he  modified  circular  slult  regisier  can  also  increase  efficiency  by  simultaneously  placing  two  or 
more  sets  ol  complex  pole  pairs  at  their  desired  location  rather  than  requiring  a  separate  circular  shift 
legister  lor  each  complex  pole  pair,  for  example  the  addition  of  a  portion  of  2N  =  22  and  2.\  - 
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6.  PERFORMANCE  CHARACTERISTICS 

6.1  LINEARITY 

Since  each  of  these  functional  building  blocks  (except  for  the  equilibration  structure)  is  structurally 
similar  to  a  stage  of  a  CCD  delay  line,  each  of  them  meets  the  requirements  of  high  speed.  Further¬ 
more,  since  each  performs  its  function  in  the  charge  domain  with  coefficients  that  depend  only  on  plan- 
view  geometry,  they  all  provide  high  linearity  and  wide  dynamic  range  regardless  of  the  choice  of  sur¬ 
face  or  buried  channel  technology.  The  only  requirement  for  accurate  operation  is  efficient,  unidirec¬ 
tional  charge  transfer.  The  speed  is  limited  only  by  the  onset  of  efficiency  loss  due  to  insufficient  time 
lor  complete  transfer  to  take  place.  This  is  a  well  studied  problem  in  the  context  of  conventional  CCD 
delay  lines,  and  may  occur  at  frequencies  of  a  few  megahertz,  for  surface  channel  structures  with  long 
gate  lengths  to  several  hundred  megahertz  for  buried  channel  devices.  Since  surface  channel  devices 
utilize  minority  carriers  and  buried  channel  devices  use  “majority  carriers”  in  a  depleted  region,  both 
of  which  are  strictly  conserved  in  the  transfer  process,  linearity  is  excellent  in  the  summation  operation 
for  a  wide  variety  of  conditions.  For  example,  even  with  buried  channel  technology,  when  two  charge 
packets  are  clocked  into  a  common  output  reservoir  the  output  charge  packet  is  equal  to  the  sum  of  the 
two  input  packets  regardless  of  their  size  or  minor  variations  in  threshold  voltage.  The  noise  sources 
that  affect  the  dynamic  range  are  kTC  noise  associated  with  the  initial  sampling  process  at  the  input  and 
transfer  noise  associated  with  incomplete  charge  transfer.  This  latter  noise  source  should  be  negligible 
for  a  high  transfer  efficiency  CCD  process.  The  secondary  effects  of  the  charge  splitting  operation,  e.g., 
splitting  accuracy,  linearity,  and  threshold  voltage  effects,  are  discussed  in  Section  6.3. 

6.2  DYNAMIC  RANGE 

The  dynamic  range  of  a  charge  domain  system  is  limited  by  the  maximum  available  potential  well 
depth  and  the  various  noise  sources.  With  the  assumption  that  kTC  noise  is  dominant,  the  signal-lo- 
noise  ratio  is  a  function  of  the  device  size,  the  clocking  voltage,  and  the  sampling  frequency.  The  max¬ 
imum  charge  packet  is 

C^mav  ~  k  ,, f  I  m.iv  ^6.  ^  ) 

and  the  noise  charge  is 

C>„„„  =  ( A  / (  ) 1  -  (6  2) 

where  I  ,,m  is  the  maximum  change  in  potential  in  a  CCD  well  before  the  charge  spills  into  adjacent 
wells,  k  is  Boltzmann's  constant  and  /  is  absolute  temperature.  For  a  typical  sized  CCD  the  input  ca¬ 
pacitance  may  be 

(  i\  -  >()()/<  m  s  I ()///;/  x  5xl()  A/>T/fj  nr  -  2.5 pi-  (6.3) 

Thus  the  input  dynamic  range  is 

s  v  =  (2,5  X  It)  '’>(2)  _ 

a,::  fl  38  x  |0  ’,)(30())(2.5  x'|()  V)\V7 

--  4  *M  x  |(F  -  <)4 ,lli  .  (6.4) 


Typical  values  of  dynamic  range  for  charge  domain  devices  assuming  no  insertion  losses  will  he 
80-100  ilH  depending  on  size.  A  larger  device  exhibits  a  wider  dynamic  range,  but  at  the  penally  of 
higher  power  consumption,  lower  device  yield,  and  higher  splitting  errors  due  to  threshold  voltage  vari¬ 
ations  depending  on  the  threshold  voltage  autocorrelation  distance. 

Another  noise  sources  to  be  taken  into  account  is  charge  transfer  noise 


t\  2  I//.Y  A  /  In? 


(6  s) 


wlicie  V  is  the  surface  stall'  density,  I  is  the  CCI)  area  pel  stage  and  ii  is  the  mmiliei  ol  liansleis 


6-1 


Sii.„r  i he  miiiilvi  til  u.insK  is  is  t\ pie. ilk  vei>  sin. ill  lui  charge  vloniam  devices.  ihc  elleet  ol  ihis  noise 
smtlu'  is  mii. i  II  I  Ik-  /■ /(  imiM  i.l  Mk-  i.ul|>iil  c liargc-lo-vullage  conversion  circuits  typically  .1  livl- 

.uni  l)u. 11  ilill umim)  v.m  Ik  taken  uul  h\  cor iclatcd  double  sampling  after  amplilicalinn  I  lie  Johnson 

in  use  11I  the  mil  pul  ampiilicr  can  he  e.ileul.ilcil  .uul  can  he  kepi  small  due  In  1  lie  small  output  capaci- 

1  a  nee  ul  1  he  ((  I)  I  he  amplifier  I/I  noise  is  also  ellcclively  cancelled  oui  hy  the  correlaled  double 

sampling  operation  when  these  devices  are  operated  al  Iheir  typical  clock  I reguencies  ol  >  I  Mil/ 

6.3  I’OKKKIC’IKNT  ACCURACY 

Minority  carrier  charge-packet  splitting  is  used  in  charge  domain  devices  to  implement  fractional 
multiplication  hy  filler  coefficients.  Since  high  precision  coefficients  are  needed  for  the  proper  imple- 
meniaiion  of  various  tillers,  ii  is  important  to  determine  and  accurately  control  1  he  factors  that 
influence  the  accuracy  and  noise  of  this  operation.' 101 

I  he  concept  of  charge-packet  splitting  is  to  divide  the  charge  in  a  CCD  register  into  portions  of  the 
total  input  charge  packet,  where  each  portion  is  manipulated  independently.  A  simple  way  for  perform¬ 
ing  this  operation  is  to  divide  the  main  CCD  channel  into  sub-channels  by  inserting  field  oxide  barriers 
into  the  channel  at  specific  locations.  The  resulting  charge-packets  are  proportional  to  the  widths  of  the 
channels  determined  hy  the  harriers. 

I  hese  harriers  can  he  constructed  in  several  different  configurations,  f  irst,  the  shape  of  the  leading 
edge  ol  the  harrier  may  he  designed  to  he  either  truncated  (blunt)  or  gradual  (point  edge).  In  the  case 
ol  the  blunt  barrier,  the  minority  charge  carriers  confronting  the  wide  edge  of  the  barrier  are  subject  to 
uncertainty  on  the  path  they  have  to  follow.  Since  there  is  a  region  where  the  carriers  have  no  clear 
choice,  and  small  uncontrolled  potential  gradients  can  become  important,  this  structure  does  not  yield  a 
high  splitting  accuracy.  Hie  structure  with  the  point-edge  barrier  has  the  CCD  channel  pinched  in  by 
half  1  he  barrier  width  from  each  side  of  the  channel  as  shown  in  Figure  6.1.  This  establishes  a  more 
symmctriial  configuration  than  the  blunt  edge  structure,  and  the  region  where  charges  can  be 
mllucnccd  hv  uncontrolled  effects  is  smaller. 

I  he  existence  ol  transfer  gates  and  storage  gates  in  a  two  phase  CCD  introduces  two  possibilities  for 
the  location  ol  1  he  barriers  II  the  leading  edge  of  the  barrier  is  located  under  the  transfer  gate  (Fig- 
ute  (>  I  »  the  charge  packet  is  split  "011-lhe-lly"  or  dynamically,  as  it  is  being  transferred  from  one 
siot.igc  gate  to  the  next  one  Since  the  potential  of  the  receiving  reservoir  is  always  larger  than  the 
tiaiisler  gate  region,  the  charge  cannot  "backup”  to  go  around  the  barriers  from  one  channel  to  the 
other.  In  this  configuration  the  split  ratio  is  determined  only  by  the  ratio  of  widths  of  the  receiving 
channels  which  is.  m  turn,  controlled  hy  the  position  of  the  point  edge  of  the  barrier. 

I  he  alternate  techim|ue  called  static  splitting  places  the  leading  edge  of  the  barrier  under  the  storage 
gale  (Figure  <>  2)  When  the  charge  carriers  arc  transferred  into  that  reservoir  they  are  distributed  uni¬ 
formly  act oss  the  entire  area,  prior  to  being  split.  When  the  transfer  gate  turns  on,  carriers  which  tire 
located  to  the  right  of  the  barrier  (to  the  right  of  line  B-  B'  in  Figure  6.2)  will  be  transferred  to  their 
roi responding  isolated  channels  I  he  initial  charge  distribution,  and  therefore  the  splitting  ratio  is 
determined  hy  the  ratio  ol  areas  above  and  below  the  barrier  region  of  Figure  6.2.  Carriers  initially  to 
the  left  ol  the  barrier,  however,  will  be  split  dynamically  and  the  split  ratio  is  therefore  determined  by  a 
iaiio  ol  channel  widths,  l  ints,  both  width  and  area  ratios  must  be  controlled  when  utilizing  the  static 
splitting  technique 

I11  order  to  obtain  a  completely  symmetrical  electric  field  pattern  at  the  gale  where  the  split  takes 
place,  a  set  of  "dummy  splitters”  is  used,  along  wilh  the  actual  coefficient  barrier.  These  splitters  are 
placed  with  identical  spacings  across  the  channel,  and  tire  only  a  few  stages  long  (Figure  6.3).  The 
numhet  ol  these  barriers  should  be  chosen  appropriately  so  that  the  incoming  charge  packet  is  divided 
into  equal  components  by  the  set  of  barriers.  Only  one  of  these  components  eventually  becomes  the 
split -oil  channel,  while  till  others  tire  recombined  after  the  few  stages.  By  creating  a  symmetrical  struc¬ 
ture.  these  dummy  splitters  serve  to  ratio  till  aspects  of  the  splitters  in  approximate  proportion  to  the 
area  I  hey  compensate  non  linearities  due  to  edge  effects  and  equalize  lateral  forces  applied  to  the 
charge  earners.  I  his  technique  is  used  both  lor  dynamic  and  static  splitting  operations. 
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F'igure  6.1.  Dynamic  barrier  charge  splitting. 

Another  method  of  achieving  charge  packet  splitting  in  CCD’s  is  to  conned  the  output  of  a  CCD 
register  to  the  input  of  two  other  CCD  registers  as  shown  in  Figure  6.4.  In  this  structure,  the  charge 
carriers  are  first  equilibrated  among  the  diodes  during  one  clock  cycle,  and  are  transferred  to  the  poten¬ 
tial  wells  of  the  split-off  channels  on  the  next  cycle.  In  this  case,  the  capacitance  ratio  of  the  input 
diffusions  determines  the  split  ratio. 

Experimental  circuits  have  been  designed  and  fabricated  to  investigate  the  packet  splitting  accuracies 
for  the  techniques  explained  above.  The  experiments  are  collected  into  a  four-chip  set  which  has  been 
fabricated  using  double  poly-silicon  gale  technology,  12.5  fim  gate  length  width  7.5  /i/ti  separation. 

To  investigate  the  dependence  of  split  accuracy  on  channel  width  and  on  coefficient  value,  split  ra¬ 
tios  designed  for  0.10,  0.25,  0.50,  0.75,  and  0.90  are  implemented  in  two  channel  widths,  250  fiw  and 
500/u/h.  The  circuits  realizing  the  0.10  and  0.90  coefficients  are  also  implemented  with  dummy 
splitters.  For  the  dynamic  splitting  technique  these  experiments  are  organized  into  two  chips,  A  and  B. 
The  designs  are  replicated  onto  chips  C  and  I)  for  the  static  splitting  experiments. 

The  “fill  and  spill”  method,  utilizing  an  input  diode  with  polysilicon  gates  for  DC  isolation,  refer¬ 
ence  and  signal  input,  is  used  for  charge  injection  into  the  registers.  The  output  of  the  CCD  consists  of 
a  precharge-and-float  circuit  and  a  source  follower,  (lie  same  output  circuit  is  used  for  the  two  (split) 
charge  packets  which  are  obtained  from  the  same  original  packet,  to  eliminate  possible  output  non- 
lineariiics  This  is  accomplished  by  having  the  register  with  the  narrow  channel  six  cells  shorter  than 
l he  ,ul|.i>  t-nl  i  li.moel  thus  pint  idmg  a  lime  delay  between  the  outputs  o!  the  two  registers 
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figure  (i.3.  Charge  splitting  using  dummy  splitter. 
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Figure  6.4.  Device  employing  input  diffusions  lor  charge  splitting. 

Tables  6.5  and  6.6  list  the  error  in  the  coefficients  with  respect  to  (he  measured  geometries  at 
1.6  MHz,  for  the  dynamic  (e  Sl/)  and  static  (e,,)  splitting  techniques,  respectively.  In  general  it  is  seen 
that  large  coefficients  are  realized  with  higher  accuracies  than  the  small  ones.  This  is  due  lo  the  fact 
that  the  small  coefficient  error  is  very  sensitive  lo  the  presence  or  absence  of  only  a  few  additional 
charge  carriers.  This  is  particularly  important  in  narrow  channel  devices,  where  the  size  of  the  charge 
packet  is  small.  The  result  is  that  split  accuracy  is  higher  in  wider  channel  devices,  especially  when 
small  coefficients  are  implemented,  l  or  example,  the  coefficient  0.10  realized  in  250  /a  m  and  500  ^  w 
channel  widths  have  e  ,,/s  of  1.76%  and  2.48%  respectively,  whereas  the  coefficient  0  90  again  reali/ed 
in  the  same  dimensions  have  e  ,,/s  of  0.29%  and  0.19%  respectively. 

The  dynamic  charge  splitting  technique  results  in  a  higher  coefficient  accuracy  than  the  static  split¬ 
ting,  due  to  the  fact  that  in  the  former  the  coefficient  is  determined  solely  by  the  barrier  location  while 
in  the  latter  it  is  also  a  function  of  an  area  ratio.  In  the  case  of  0.25  and  0.75  coefficients  with  500  fim 
channel  widths  and  dynamic  splitting,  the  ew’s  are  0.41%  and  0.13%,  while  the  0.2459  and  0.7541 
coefficients  which  use  static  splitting  have  « ,,/s  of  0.79%  and  0.25%  respectively. 

The  use  of  dummy  splitters  significantly  increases  the  coefficient  accuracy.  For  example  in  the  case 
of  1:9  dynamic  split,  implemented  with  500  ^ m  channel,  the  introduction  of  dummies  reduces  the 
effect  from  1.76%  to  0.73%  for  the  0.10  coefficient  and  from  0.19%  to  0.07%  for  the  0.90  coefficient. 
Blunt  splitters,  however,  decrease  the  measured  accuracy.  Again  for  a  500 /am  channel  and  a  1:3  static 
split,  the  coefficient  error  increases  by  about  a  factor  of  three  when  the  blunt  rather  than  sharp  barrier 
is  used.  Table  6.7  lists  the  coefficient  accuracies  for  dynamic  splitting  and  500 /a  m  channel  at  1.6  and 
5.4  MHz.  Comparison  ol  these  values  with  those  at  73.5  KHz  (  Table  6.8)  shows  that  errors  are  almost 
consistent  up  to  1.6  MHz  and  they  slightly  increase  at  5  4  Mllz. 

The  devices  utilizing  input  diffusions  for  charge  splitting  do  not  provide  as  high  accuracy  as  the  oth¬ 
ers.  Since  the  charge  is  split  during  the  on-time  of  the  gates  next  to  the  receiving  diodes,  the  field  ox¬ 
ide  separating  these  two  diffusions  acts  as  a  blunt  splitter.  Tor  a  500  /am  channel  the  0.25  and  0.75 


Table  6.5 

SUMMARY  OK  THE  RESULTS  FOR  DYNAMIC  SPLITTING 


Split 

Ratio 

Design 
Channel 
Width  (/cm) 

Design 

Coeff. 

Measured 

Geometry 

Measured 

Signal 

«  w 

1:1 

125.0 

0.50 

0.5020 

0.5042 

0.44% 

125.0 

0.50 

0.4980 

0.4958 

0.44% 

62.5 

0.25 

0.2545 

0.2481 

2.52% 

1:3 

187.5 

0.75 

0.7455 

0.7419 

0.86% 

w/o 

D's 

125.0 

0.25 

0.2443 

0.2453 

0.41% 

375.0 

0.75 

0.7557 

0.7547 

0.13% 

25.0 

0.10 

0.0926 

0.0903 

2.48% 

1:9 

225.0 

0.90 

0.9074 

0.9100 

0.29% 

w/o 

Ds 

50.0 

0.10 

0.0964 

0.0947 

1.76% 

450.0 

0.90 

0.9036 

0.9053 

0.19% 

25.0 

0.10 

0.0926 

0.0938 

1.30% 

1:9 

225.0 

0.90 

0.9074 

0.9062 

0.13% 

w.D’s 

50.0 

0.10 

0.0964 

0.0957 

0.73% 

450.0 

0.90 

0.9036 

0.9042 

0.07% 

62.5 

0.25 

0.2545 

0.2255 

1 1.59% 

1:3 

187.5 

0.75 

0.7455 

0.7745 

3.89% 

w/ln. 

Diff. 

125.0 

0.25 

0.2443 

0.2401 

1.72% 

374.0 

0.75 

0.7557 

0.7599 

0.56% 

25.0 

0.10 

0.0926 

0.0665 

28.19% 

1:9 

225.0 

0.90 

0.9074 

0.9335 

2.88% 

w/ln. 

Diff. 

50.0 

0.10 

0.0964 

0.0802 

16.80% 

450.0 

0.90 

0.9036 

0.9198 

1.79% 

coefficients  are  realized  in  ihis  ease  wilh  1.72%  and  0.56%  errors  which  are  about  2  limes  higher  than 
the  corresponding  values  in  the  ease  of  static  splitting  and  a  blunt  barrier. 

The  coefficient  accuracies  obtained  in  the  experiments  are  generally  high  enough  to  build  high  per¬ 
formance  signal  processing  devices,  and  the  major  error  sources  can  in  some  cases  be  compensated. 
C  ompurison  of  the  design  and  physical  geometries  indicate  that  a  substantial  amount  of  error  is  intro¬ 
duced  In  process  bias.  Since  physical  ge*. metrics  are  determined  by  the  field  oxide  layer,  the  lateral  ox¬ 
ide  growth  and  ihc  isotropic  etching  of  the  nitride  layer  used  to  define  the  field  oxide  can  introduce  er¬ 
rors  into  the  geometries.  I  he  former  can  be  compensated  by  calibration  of  this  effect  and  predislorting 
the  field  oxide  mask.  The  latter  can  be  minimized  by  using  anisotropic  etching  techniques  to  minimize 
the  intrinsic  potential  well  nonuniformiiy  at  the  edges  of  the  channel.  The  so-called  edge  effect,  is  the 
second  error  source.  The  dummy  splitters  used  to  minimize  this  effect  appear  to  have  made  a 
significant  improvement,  frequency  limitation  of  the  drivers  is  the  present  limitation  in  very  high  fre¬ 
quency  operation  This  in  turn  can  be  substantially  improved  with  more  suitable  components  and  Rf 
design  techniques. 


SUMMARY  OF  THE  RESULTS  FOR  STATIC  SPLITTING 


Design 


Split 

Ratio 

Channel 
Width  (nm) 

Design 

Coeff. 

Measured 

Geometry 

Measured 

Signal 

€  */ 

1:1 

125.0 

0.5000 

0.5020 

0.4954 

1.31% 

125.0 

0.5000 

0.4980 

1.33% 

62.5 

0.2417 

0.2464 

0.2330 

5.44% 

1:3 

187.5 

0.7583 

0.7536 

0.7670 

1.78% 

w/o 

D’s 

125.0 

0.2459 

0.2405 

0.79% 

375.0 

0.7541 

0.7595 

0.7614 

25.0 

0.0868 

0.0792 

1:9 

225.0 

0.9132 

w/o 

D’s 

50.0 

0.09356 

0.0898 

4.23% 

mmm 

25.0 

BBSS 

0.0927 

5.18% 

1:9 

225.0 

w/D’s 

50.0 

0.0964 

1.87% 

1:3 

125.0 

0.2415 

0.2469 

2.24% 

Blunt 

375.0 

0.7585 

Table  6.7 

COEFFICIENT  ERRORS  VS.  FREQUENCY 


Coe  IT. 

€.v</  at 

1.6  MHz 

e  v,/  at 
5.4  MHZ 

0.10 

1.76% 

2.70% 

Without 

0.25 

0.41% 

0.98% 

Dummies 

0.75 

0.13% 

1.32% 

0.90 

0. 1 9% 

0.29% 

With 

0.10 

0.73% 

1.87% 

Dummies 

0.90 

0.07% 

0.20% 

Table  6.8 

COEFFICIENT  ERRORS  AT  73.5  KHz 


Design 


Split 

Ratio 

Channel 
Width  ifim) 

Design 

Geometry 

Measured 

Geometry 

Measured 

Signal 

€*/ 

1:1 

125 

0.50 

0.5020 

0.5029 

0.18% 

125 

0.50 

0.4980 

0.4971 

0.18% 

1:3 

62.5 

0.25 

0.2545 

0.2429 

4.56% 

w/o 

187.5 

0.75 

0.7455 

0.7571 

1.56% 

D’s 

1:3 

125 

0.25 

0.2443 

0.2412 

1.27% 

w/o 

375 

0.75 

0.7557 

0.7588 

0.41% 

D's 

1:9 

25 

0.10 

0.0926 

0.0877 

5.39% 

w/o 

D's 

225 

0.90 

0.9074 

0.9123 

0.55% 

1:9 

50 

0.10 

0.0964 

0.0933 

3.22% 

w/o 

D’s 

450 

0.90 

0.9036 

0.9067 

0.34% 

1:9 

25 

0.10 

0.0926 

0.0919 

0.86% 

w/D’s 

225 

0.90 

0.9074 

0.9081 

0.08% 

1:9 

50 

0.10 

0.0964 

0.0970 

0.62% 

w/D’s 

450 

0.90 

0.9036 

0.9030 

0.97% 
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7.  TESTING  OF  CHARGE  DOMAIN  DEVICES 

The  testing  and  characterization  of  charge  domain  devices  can  he  performed  cither  by  looking  at 
time  domain  or  frequency  domain  analyses.  Time  domain  is  the  chosen  technique  for  this  work  and  a 
computerized  data  acquisition  system  has  been  developed  to  analy/e  the  results.  The  computer  analysis 
program  described  in  the  next  section  is  used  to  relate  the  time  domain  analysis  to  frequency  domain 
results 

7.1  TESTING  PHILOSOPHY 


Although  the  ultimate  goal  in  designing  charge  domain  filters  is  to  reali/e  a  certain  frequence 
response,  the  devices  themselves  are  constructed  as  structures  which  manipulate  charge  packet  embodi¬ 
ments  of  samples  of  analog  data  in  the  time  domain.  Thus,  the  expected  performance,  and  any  devia¬ 
tions  therefrom  can  best  be  understood  by  looking  at  the  functioning  of  the  charge  domain  building 
blocks  in  the  time  domain  when  possible.  In  this  framework  inaccuracies  in  the  impulse  response  van 
be  directly  related  to  performance  due  to  functions  such  as  charge  splitting  errors,  charge  transfer 
inefficiency  effects,  or  errors  in  the  initial  design.  This  technique  has  given  a  deeper  understanding  of 
some  of  the  error  sources  in  certain  designs  and  has  led  to  the  development  and  demonstration  of  im¬ 
proved  designs.  It  must  always  be  understood  however,  that  the  final  purpose  is  usually  frequency 
based,  and  as  such,  the  importance  of  time  domain  errors  or  noise  must  be  related  to  frequency  domain 
performance.  This  transformation  can  he  made  using  computer  simulation  tools  and  is  described  later. 

Another  advantage  of  using  the  impulse  or  pulse  response  data  is  that  averaging  over  many  cycles  of 
the  data  may  be  done  to  improve  the  accuracy  by  averaging  out  the  noise  effects.  The  averaging  opera¬ 
tion  is  used  to  test  a  wider  dynamic  range  than  may  be  available  with  a  single  sample  through  an 
analog-to-digital  (A/D)  converter.  If  the  device  noise  and  the  A/I)  quantization  noise  are  uncorrelated 
from  sample  to  sample,  and  the  signal  is  fully  correlated,  i.e.,  no  drift  terms,  then  the  signal  to  noise 
improves  with  the  square  root  of  the  number  of  readings.  With  this  technique  terms  such  as  splitting 
accuracy  can  be  measured  to  very  high  precision.  Meanwhile  by  taking  the  standard  deviation  of  these 
readings,  the  various  noise  sources  themselves  can  be  analyzed. 

7.2  TEST  STATION  DESIGN 

A  test  station  has  been  designed  and  built  which  allowed  characterization  of  the  charge  domain  de¬ 
vices.  This  station  consisted  of  two  parts,  the  electronics  to  run  the  chip  itself  and  a  microcomputer 
with  data  acquisition  capabilities  to  acquire  and  analyze  the  data.  A  block  diagram  for  the  driver  elec¬ 
tronics  is  shown  in  figure  7.1.  It  consists  of  a  voltage  controlled  oscillator  for  variable  clock  speeds, 
some  TTL  circuits  to  generate  the  proper  timing  signals,  level  shifters  to  generate  the  MOS  clocking 
levels  needed  for  the  chip,  and  variable  DC  supplies  for  the  chip.  In  addition  there  is  a  differential 
amplifier  on  the  output  and  a  sample  and  hold  device.  The  sampled  and  held  output  is  available  for  ac¬ 
quisition  by  the  computer. 

The  computer  includes  a  Z-80  central  processing  unit,  an  RS-232  terminal  interface,  a  Basic  inter¬ 
preter  in  ROM  with  additional  ROM  for  the  Basic  analysis  programs  and  the  operating  system  software. 
The  computer  also  contains  program  RAM  plus  IhK  of  RAM  storage  for  the  data.  A  CRT  or  printer 
terminal  can  be  connected  for  input/oulput.  A  custom  A/D  conversion  board  has  been  built  for  the 
system.  This  board  allowed  acquisition  of  the  data  with  12  bits  of  accuracy  at  rates  up  to  I  Mllz.  This 
hoard  included  the  necessary  timing  for  loading  the  data  directly  into  the  RAM  storage.  Once  loaded, 
the  data  was  available  for  analysis  by  the  Basic  analysis  programs.  One  such  program  is  listed  in  Appen¬ 
dix  I.  T  his  program  is  used  to  print  out  the  time  samples  of  an  impulse  or  pulse  response  of  a  filter 
section  (see  figure  7.2).  The  decaying  sine  wave  impulse  response  can  be  noted  starting  at  data  point 
II  Averaging  of  up  to  100  readings  (#  READS  parameter)  can  be  used  to  improve  the  accuracy  of 
these  readings.  The  standard  deviation  of  these  readings  is  also  printed  out  to  give  a  measure  of  the 
noise  on  the  readings  and  also  to  assure  that  the  device  noise  is  larger  than  the  A/D  quantization  error. 
I  Ins  tool  pi ov es  to  he  quite  valuable  lor  measuring  coefficient  accinaev.  svsiem  noise,  and  transfer 
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DEVICE  TESTING 


Figure  7.1.  Test  electronics. 
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Typical  Output 


Z80-Monitor 

>  B 

#  Reads,  #  Displayed?  4,32 


Data 
Point  # 

Output 

Volts 

Noise 

#Volts 

Output 
Maxout=  1 

Noise 
Maxout  =  1 

1 

-6.10352E-04 

1.72634E-03 

—  1.11 545l-:-03 

3. 15497  E-03 

2 

—  3.051 76E-04 

1.49505E-03 

— 5.57725K-04 

2.73228E-03 

3 

-1.2207E-03 

3.07204E-03 

-2. 2  809  E-03 

5.6143E-03 

4 

0 

2.1 1432E-03 

0 

3.8640E-03 

5 

-9. 15527E-04 

3.62804E-03 

— 1.673 17E-03 

6.63044E-03 

6 

— 9.1 5527E-04 

1.79683E-03 

— 1.673 17E-03 

3.28379E-03 

7 

— 3 .05 1 76E-24 

2.49175E-03 

—  5.57725E-04 

4.5538E-03 

8 

—  1 .2207E-03 

2.90586E-03 

—  2.2309E-03 

5.3106E-03 

9 

0 

2.1 1432E-03 

0 

3. 86403  E-03 

10 

— 3.051 76E-04 

2.28373E-03 

-5.57725E-04 

4.17363E-03 

11 

.115356 

2.63702E-03 

.21082 

4.81929E-03 

12 

.279236 

0.72958E-03 

.510318 

4.98844E-03 

13 

.426941 

2.33747E-03 

.780257 

4.271 85E-03 

14 

.518188 

9.967E-04 

.947016 

1 .82 1 52E-03 

15 

.54718 

1.22L07E-03 

1 

2.2309E-03 

16 

.515442 

2. 1 1432E-03 

.941997 

3.86403E-03 

17 

.440674 

1.57592E-03 

.805354 

2.88008E-03 

18 

.338135 

2.90586E-03 

.617959 

5.3106E-03 

19 

.239868 

2.63702E-03 

.438372 

4.8I929E-03 

20 

.146179 

3.48845E-03 

.26715 

6.37532E-03 

21 

.0726318 

1.9934E-03 

.132738 

3.64304E-03 

22 

.0177002 

1. 40955  E-03 

.032348 

2.57602E-03 

23 

-.0146484 

2.541  IE-03 

-.0267708 

4.64398E-03 

24 

-.0292969 

2.90586E-03 

-.0535416 

5.3106E-03 

25 

-.0299072 

1.72634E-03 

-0.D54657 

3.15497E-03 

26 

-.0213623 

3.30568E-03 

-.0390407 

6.041 3  E-03 

27 

-.0128174 

1.72634E-03 

-.0234244 

3.15497E-03 

28 

— 9. 15527E-04 

3.48845E-03 

— 1.673 17E-03 

6.37532E-03 

29 

7.01904E-03 

2.28373E-03 

.128277 

4. 17363E-03 

30 

.012207 

2.1 1432E-03 

.022309 

3.86403E-03 

31 

.015564 

3.34303E-03 

.028444 

6. 19057E-03 

32 

.0152588 

9.967E-04 

.0278862 

1.82I52E-03 

RETAKE  DATA  (Y  OR  N)?  N 
REPEAT  CALCULATION  (Y  OR  N)?  N 
OK 


Figure  7.2.  Measured  impulse  response  of  charge  domain  low  pass  filler. 


8.  COMPl .TER  SIMULATIONS  TOOLS 


Once  the  time  domain  data  is  taken,  it  is  useful  to  compare  the  experimental  data  with  the  design 
data.  When  the  application  involves  tillering,  the  specifications  are  usually  given  in  terms  of  frequency 
response.  Therefore  computer  tools  are  required  to  convert  ihc  lime  domain  data  into  die  frequency 
domain  so  that  the  experimental  frequency  response  may  he  determined  and  compared  with  the  desired 
response.  Computer  tools  are  also  required  w  hen  undertaking  i he  design  ol  a  new  tiller  I  he  proper 
placement  of  poles  and  /e roes  within  the  constraints  imposed  In  charge  domain  technology  can  he  ac¬ 
complished  through  simulation,  from  these  simulations  it  can  he  determined  whether  a  lilici  design  is 
practical  with  charge  domain  technology. 

The  computer  program  "DSP"  has  been  written  lo  address  these  needs  \  listing  is  given  in  Ap¬ 
pendix  2.  It  consists  of  a  number  ol  LOR  I  RAN  (>(>  suhiouiiues  written  on  a  Honeywell  <>1)00  main¬ 
frame  computer.  The  routines  cover  functions  used  for  digital  (or  discrete  time*  signal  processing  as 
well  as  printing  and  plotting  routines  for  visualizing  the  outputs.  The  signal  processing  routines  include 
routines  like  ZPOLR  which  finds  the  roots  of  a  polynomial,  P/.IMP  which  calculates  the  impulse 
response  from  the  pole-zero  locations,  LRLQI*/  which  converts  zero  locations  stated  as  fractions  ol  the 
clocking  frequency  to  /.-plane  zero  locations,  PX  which  calculates  the  frequency  response  from  the  pole 
and  zero  locations,  and  subroutine  COLT  which  calculates  the  frequency  response  from  die  numerator 
and  denominator  coefficients  of  a  transfer  function.  Other  routines  include  LIT  which  performs  a  last 
Lourier  transform  for  converting  between  time  and  frequency  domain  results,  MAO  which  converts  a 
complex  frequency  response  into  magnitude  (linear  and  dB>  and  phase  information,  and  SAMPLL 
which  multiplies  the  frequency  response  data  bv  the  si n ( x ) / x  response  of  a  data  sampler.  Printing  rou¬ 
tines  include  PRNTIMP,  PRNTLRLQ,  and  PRNTPZ  for  impulse  response,  frequency  response  and 
pole-zero  information  respectively,  as  well  as  P/MINMAX  which  prints  out  the  local  minimums  and 
maximums  of  frequency  response  curves.  The  plotting  routines  are  PI  LIMP,  Pl.TMAO.  PI  TPIIASL. 
and  PLTP/  for  impulse  response,  magnitude  and  phase  and  pole-zero  plots.  The  subroutine  PI  T  is  the 
general  purpose  X-Y  plotting  routine  called  by  Pl.TIMP,  Pl.TMACi,  and  PI. TPIIASL.  These  routines 
call  subroutines  of  Hewlett-Packard's  PLOT21  plotting  package  which  interlaces  with  I  IP's  722 1  A  pen 
plotter. 

With  the  availability  of  these  routines,  a  simple  program,  such  as  the  one  listed  in  Appendix  3  and 
excerpted  in  f  igure  8.1,  is  all  that  is  required  to  design  a  charge  domain  filler.  The  listed  program  was 
used  to  design  the  8  pole  bandpass  filler  which  has  been  designed  under  this  contract.  There  is  one 
section  for  each  of  the  four  pairs  of  complex  conjugate  poles,  and  one  section  for  the  cancellation 
zeroes.  Each  section  may  be  analyzed  by  itself  or  in  combination  with  any  or  all  of  the  various  other 
sections.  Lor  each  of  the  pole  sections  the  polynomial  coefficients  must  be  specified.  These  correspond 
directly  to  the  splitting  coefficients  where  the  polynomial  subscript  (actually  the  exponent  of  /  '  I 
corresponds  to  the  number  of  stages  in  the  feedback  loop.  Lor  the  bandpass  Tiller  design  described  in 
Section  9.  the  0  =  91"  poles  use  a  charge  domain  stage  which  implements  the  function 

//(/)  =  .1/(1  +  .(>/  14  +  3Z  l7)  <5.1 ) 

The  required  polynomial  coefficients  for  the  program  tire  PCOLLd)  -  I,  IK'OLIT 1 5)  —  h.  and 
PC()LL(I8)  =  .3.  The  corresponding  charge  domain  stage  has  a  feedback  portion  =  h  with  A'  =  14 
CCD  stages  of  delay,  another  feedback  portion  =  .3  with  .V  =  17  stages,  and  the  remaining  portion  = 
.1  sent  directly  lo  the  stage  output.  The  DSP  program  produces  pole-zero,  frequency  response,  and 
time  domain  plots,  as  well  as  printouts  of  pole  and  zero  locations  and  filler  coefficients.  Examples  of 
these  plots  and  printouts  tire  displayed  throughout  Section  9.  The  feedback  coefficients  and  the  number 
of  stages  in  the  feedback  loop  may  be  varied  and  the  effect  on  the  response  can  he  visualized  in  the 
pole-zero  and  frequency  response  plots.  Thus,  each  of  the  pole  sections  can  be  optimized  under  the 
charge  domain  constraints. 

Once  the  remaining  poles  have  been  placed  to  shape  the  passband,  the  cancellation  zeroes  can  he 
pl.uvd  to  cancel  out  the  spurious  poles  Ibis  is  done  bv  placing  the  transversal  filler  zeroes  on  the  unit 
>iul<  ii  iId.-  1 1 1  spirm  ics  ol  iId-  spmioiis  poles  I  Ins  is  deiuonsl  I  tiled  ill  the  lisling  ol  I  igiue  8  I  vcltcie 
ill)-  zeio  l*ii  .i  1 1)  >  i  is  ,iic  spe>  tiled  .is  positions  on  the  unit  cilcle  as  a  li.ulion  ol  the  distance  aiound  I  lu¬ 
ll  1 1 1 1  encle  <ic  as  a  I  ruction  ol  the  clocking  frequency),  finally,  the  oveiall  frequency  response  and 
pole  /cm  plots  ol  the  filter  design  can  he  displaced  and  compared  with  the  filler  specifications. 


10  print  *,  ’ 

print  \  '91  DEG  POLES' 
NPOLE=  17 
DO  12  1  =  l,NPOLE+  1 
12  PCOEF(l)  =0. 

PCOEF(l)  =  1. 
PCOEF(15)  =  .6 
PCOEF(18)  =  .3 
KLAB  =  30 
CALL  CALCPOLE 
DO  14  1  =  1, N 
14  GO)  =G(I)*H(I) 


C  40  print  ’ 

print  *, 'ZEROS’ 

NPOLE=0 
NFREQ  =  6 

DATA  (ZFREQ(I),I  =  1,6)  /.0325,. 108,. 180,. 320,. 392, .4675/ 
KLAB  =  20 
CALL  CAL  CZERO 
DO  46  1  =  1  .N 
46  GO)  =G(1)*H(I) 

C 

^  ^  i****************:*#**********************1*1** 


Figure  8. 1  Partial  listing  of  bandpass  filter  design  program. 

I  his  program  has  been  invaluable  in  the  development  of  design  methodologies  lor  charge  domain 
lev  ices,  and  in  analysis  ol  their  experimental  performance.  The  program  has  recently  been  upgraded  to 
OR  I  RAN  77  on  a  Digital  Equipment  Corporation  VAX  11/782  computer.  The  plotting  routines  have 
ilso  been  enhanced  with  such  Features  as  windowing  and  curve  smoothing  and  allow  plotting  on  Tek- 
ronix  graphics  terminals  as  well  as  Hewlett-Packard  pen  plotters.  Thus,  quick  visual  Feedback  is  uvail- 
ihle  on  a  CRT  terminal,  with  hard  copy  plotting  once  the  design  has  been  confirmed. 
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1  lie  second  goal  of  the  contract  as  listed  in  Section  I,  was  the  design  and  fabrication  of  a  test  vehi¬ 
cle  tiller  which  demonstrates  the  unique  capabilities  of  the  charge  domain  approach  Before  embarking 
on  the  design  ol  a  lest  vehicle  filter,  ;m  analysis  of  an  existing  charge  domain  low  pass  liltci  i  hip 
"(  Oi  l  was  undertaken.  I  hese  results  led  to  a  decision  that  a  redesign  ol  the  low  pass  chaigc  domain 
filter  would  provide  further  insight  into  design  considerations.  This  redesigned  chip  entitled  (l)l.1.  ex 
h  i  hi  ted  much  improved  performance  over  our  earlier  design,  following  this  accomplishment,  the 
design  of  a  high  performance  bandpass  filter  was  attempted  This  exercise  led  to  .1  thorough  undei- 
standmg  of  the  design  tradeoffs  which  optimize  charge  domain  filter  designs  An  8-polc  narrow  hand 
filter  was  chosen  as  the  test  vehicle,  because  its  simulated  performance  was  unachievable  with  alterna¬ 
tive  technologies.  The  design  was  implemented  in  silicon  and  the  resulting  chips  entitled  (  DM 
demonstrated  a  frequency  response  which  matched  the  simulations  very  closely. 

9.1  CHARGK  domain  low  pass  fii.ikr 

9.1.1  First  results  -  C'DFl 

I  he  design  concepts  described  in  Section  2  have  been  utilized  in  the  design  of  the  first  three  pole 
charge  domain  filter  "('Dl  l.''  A  Butterworlh  low  pass  filter  has  been  chosen  for  simplicity  of  design 
while  still  demonstrating  the  charge  domain  concepts  of  charge  splitting,  equilibration,  feedback  and 
cascading  of  sections.  The  design  requires  a  recursive  accumulator  to  implement  the  real  axis  pole  and 
the  complex  pile  pair  is  implemented  by  a  circular  shift  register  in  combination  with  a  cascaded 
transversal  filter  to  provide  the  cancellation  zeroes.  The  chip  has  been  layed  out  and  fabricated  in  the 
(if  processing  facility  in  Schenectady,  NY. 

A  photomicrograph  of  the  chip  is  shown  in  f  igure  9.1.  The  chip  contains  two  fill-and-spill  input 
sections  (positive  and  negative  inputs),  two  recursive  accumulators  to  implement  the  real  pole,  a 
16-stage  recursive  loop  to  provide  the  8  poles  around  the  unit  circle,  two  7-slage  transversal  filters  to 
cancel  out  6  of  those  8  poles,  and  two  output  sections  (positive  and  negative  output).  The  device  litis 
been  processed  in  7.5  fxm  p-channel  technology  (12.5  /cm  long  storage  gates  with  7.5  /uni  long  transfer 
gates'.  The  chip  size  is  55  x  88  mil’  including  the  large  area  used  for  bonding  pads.  Chip  power  is  less 
than  10  mW  at  a  I  Mllz  clock  frequency.  A  pseudo-one  phase  clocking  scheme  is  used  to  move  the 
charge  packets  and  10  n  wide  thick  oxide  barriers  (narrowed  to  a  point)  are  used  to  split  the  charge 
packets. 

A  number  of  test  devices  are  included  in  the  mask  set  so  that  each  of  the  subsections  of  the  tiller 
can  be  tested  individually.  The  time  response  of  several  of  these  subsections  is  shown  in  figure  9.2. 
I  igure  9.2a  shows  the  impulse  response  of  the  transversal  filter  alone.  This  is  a  finite  impulse  response 
filter  with  7  non-zero  values  as  can  be  seen  from  the  f  igure.  When  the  circular  shift  register  is  includ¬ 
ed,  the  response  continues  as  a  decaying  sine  wave  for  an  infinite  number  of  cycles  as  shown  in 
f  igure  9.2b  and  9.2c.  The  impulse  response  of  the  real  pole  section  of  f  igure  9. 2d  is  seen  to  be  the  ex¬ 
pected  decaying  exponential.  The  impulse  response  of  the  entire  filler  is  shown  in  f  igure  9.2e. 

These  time  domain  responses  have  been  analyzed  in  detail  by  the  computer  data  acquisition  system 
described  in  Section  7  and  (he  simulation  tools  described  in  Section  8.  This  analysis  enables  a  number 
of  design  problems  to  be  isolated.  While  the  device  demonstrates  the  feasibility  of  the  charge  domain 
approach,  these  problems  limits  its  performance  in  several  ways. 

A  design  error  causes  the  feedback  coefficient  for  the  recursive  loop  to  be  too  small,  thus  giving  in¬ 
complete  cancellation  of  the  excess  poles  by  the  zeroes.  liven  when  this  error  is  taken  into  account, 
the  measured  feedback  coefficient  is  in  error  and  showed  poor  linearity.  This  problem  is  caused  by  the 
equilibration  splitting  technique  itself.  Whenever  equilibration  splitting  is  to  be  used  (either  gated  or 
non-galed)  great  care  must  be  taken  to  properly  ratio  all  significant  splitting  parameters,  namely  storage 
area  oxide  capacitances,  storage  area  diffusion  capacitances,  channel  widths,  and  gale  area  itself  I  his 
lopu  is  l-ii-  desciibed  m  mine  detail  in  the  section  on  liltci  redesign 

I  he  dinipc  domain  liunxveisnl  liltci  has  been  analyzed  m  gieul  delad,  and  the  icsulls  aic  listed  m 
table  9  .!  |  hois  in  i|v  ailual  1  espouse  aic  on  the  order  ol  1%  when  compared  to  the  design  values 

However  when  the  geometrically  measurer!  coefficients  tire  compared  to  the  electrically  measured 
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f  igure  9. 1.  C'DFI  photomicrograph. 

coefficients,  the  agreement  is  much  closer.  The  discrepancy  between  the  designed  and  the  actual  chan¬ 
nel  widths  is  due  to  geometry  alterations  which  occur  during  device  processing.  This  error  has  been  re¬ 
duced  in  subsequent  devices  by  predistorting  the  mask  design  to  compensate  for  this  effect.  With  this 
compensation  and  other  improvements  described  in  Section  6.3,  the  dynamic  charge  splitting  technique 
provides  a  very  accurate  and  linear  means  for  determining  coefficients. 

I  he  configuration  of  the  real  pole  (f  igure  9.4)  section  causes  a  number  of  problems.  In  this  ap¬ 
proach  a  charge  packet  enters  storage  area  A.  The  equilibration  gale  is  turned  on  causing  the  surface 
potential  of  storage  areas  A  and  It  to  equilibrate.  When  the  gate  is  turned  off  the  charge  should  be  split 
m  the  same  ratio  as  the  areas.  Unfortunately  the  charge  under  the  equilibration  gate  does  not  split  by 
the  proper  ratio  as  the  gate  is  turned  off  but  usually  splits  closer  to  a  1:1  ratio  because  of  the  charge 
under  the  gate  region  itself,  which  lends  to  split  equally  as  the  gale  is  turned  oil'.  This  effect  causes  a 
3""  error  in  the  coefficient  of  the  real  axis  pole.  A  solution  for  this  problem  is  described  in  the  next 
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(d)  Recursive  accumulation 
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Another  problem  with  the  configuration  of  Figure  9.4  is  the  long  gates  of  the  output  section  C.  For 
:  saturation  charge  handling  capabilities  the  area  under  electrode  C  must  equal  the  area  under  electrode 

A,  and  since  the  width  is  confined  the  length  is  increased  from  12.5  to  30  u  rn.  The  transfer 

l  inefficiency  of  these  long  stages  causes  the  clocking  speed  of  these  devices  to  be  limited  to  less  than 

1  MHz  (Figure  9.5).  This  problem  has  also  been  overcome  in  the  redesign.  It  is  noted  that  the  lateral 
equilibration  using  metallized  diffusions  was  very  fast  and  does  not  provide  any  speed  limitation  for 
[  these  devices. 

£  Finally  it  is  noted  that  harmonic  distortion  is  quite  low  -55  to  -60  dB  due  to  using  differential  (pos- 

t  itive  and  negative)  signal  processing  (Table  9.6). 

r 

£  9.1.2  The  Redesign  -  CDF2 

r  A  redesign  of  CDF1  has  been  completed  in  order  to  optimize  the  performance  of  the  low  pass  filter, 

1.  and  to  overcome  the  faults  which  are  apparent  in  the  operation  of  CDF1.  The  three  pole  Butterworth 

S  filter  characteristic  is  improved  by  moving  the  poles  to  the  following  locations: 


(e)  Entire  filter 
Figure  9.2.  (Cont’d) 


Table  9.6 


HARMONIC  DISTRIBUTION  FOR  A  CLOCK  FREQUENCY  OF  100  1  Hi 
AND  A  SIGNAL  FREQUENCY  OF  1  kHz 


F  Out 
(kHz) 

V  Out 
(mV) 

VOut 

(mV) 

V  Out 
(mV) 

V  Out 
(mV) 

VOut 

(mV) 

1 

10 

100 

300 

1000 

3000 

2 

.0045 

.004 

.15 

.32 

7.8 

3 

.0040 

.015 

.005 

.17 

950 

4 

.0035 

.008 

.06 

.08 

3.0 

5 

.0030 

.008 

.005 

.06 

5.1 

(Note:  all  voltages  are  rms  values) 

P,  =  0.714 

P2  =  0.778  4-  j  0.332 

Pi  =  0.778  -  j  0.332 

This  provides  an  equi-ripple  type  frequency  response  with  less  passband  attenuation,  greater  stopband 
attenuation  and  a  sharper  transition  region  (Figure  9.7). 


An  attempt  has  been  made  to  overcome  each  of  the  deficiencies  of  the  previous  niter  with  the  main 
goals  being  improved  accuracy  and  speed  performance.  This  device  has  been  produced  with  n-channel 
technology  and  with  design  rules  shrunk  by  20%.  This  yielded  a  10  yun  storage  gate  length  and  6 
transfer  gate  length.  The  saturation  charge  capability  is  therefore  reduced  and  the  low  frequency 
transfer  inefficiency  increases,  but  the  maximum  clocking  speed  should  be  increased.  The  pseudo-one 
phase  clocking  scheme  is  maintained  and  the  thick  oxide  barrier  width  is  shrunk  to  8  /am.  A  number 
of  design  changes  have  been  incorporated  and  the  improvements  have  been  successful.  The  finished 
chip  photomicrograph  is  shown  in  Figure  9.8. 

This  time  dynamic  splitting  with  barriers  is  used  to  determine  the  feedback  coefficient  for  the  recur¬ 
sive  loop,  rather  than  the  equilibration  scheme  used  for  CDF1.  This  leads  to  much  improved 
coefficient  accuracy  and  linearity.  The  old  “form  -  C”  structure  and  new  “Figure  —  8”  structure  ap¬ 
proaches  are  compared  in  Figure  9.9.  Mask  predislortion  is  used  and  coefficient  accuracy  for  the  feed¬ 
back  coefficient  (0.05%  error)  and  the  transversal  filter  coefficients  (0.1  —  0.2%)  is  much  improved 
(see  Table  9.10). 

Two  different  techniques  of  overcoming  the  deficiencies  of  the  real  pole  section  (shown  again  in 
Figure  9.1la)  have  been  evaluated.  The  first  one  is  once  again  a  gated  equilibration  technique.  Here, 
however,  the  incoming  charge  packet  is  dumped  into  the  storage  area  B  rather  than  A  (see 

CHARGE  DOMAIN  FILTERS  —  CDF2 
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Figure  9.8.  CDF2  photomicrograph. 


C"  siruciure  tor  implementing  complex  pole  pairs. 
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Table  9.10 

CI)F2  TRANSVERSAL  FILTER  COEFFIC  IENT  ACCURACY 


CoeF. 

Theory 

Actual 

Si/e 

100  kll  / 
Exp. 

29  Mil/ 
Exp. 

1 

.585 

.572 

.565 

1 

LT, 

2 

.909 

.906 

.915 

.909 

3 

1.000 

1.000 

1.000 

1.000 

4 

.911 

.910 

.916 

.909 

5 

.709 

.704 

.705 

.721 

6 

.456 

.451 

.458 

.472 

7 

.208 

.202 

.206 

.221 

FB 

1.204 

1.207 

1.200 

1.209 

EQUIL.  GATE 


30m 

(a)  Equilibration  technique  used  for  C'DFl 
10m 


L-<. 


EQUIL.  GATE 

(b)  Equilibration  technique  used  lor  CDF2 


10M 

<c>  Barrier  charge  splitting  technique  used  for  CDF2 
Figure  9. 1 1  Recursive  accumulator  structures. 
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Figure  9.11b)  This  necessitates  ;i  wider  structure  to  accomodate  the  saturation  signal.  hut  it  does  not 
incorporate  any  long  gates  and  thus  is  much  faster  than  the  earlier  approach  In  addition  great  care  lias 
been  taken  to  properly  ratio  all  capacitances  and  to  properly  ratio  the  amount  of  stored  charge  undei  the 
equilibration  gate  as  it  is  turned  oil'.  The  resulting  device  shows  high  accuracy  <0.2%  error)  and 
operates  at  high  speed  (to  >  10  Mil/). 

The  second  approach  is  an  attempt  to  use  dynamic  (barrier)  splitting  to  determine  the  coefficient  A 
barrier  split  with  a  single  stage  feedback  loop  is  implemented  as  shown  in  Figure  0.1  U\  Mere,  however, 
lateral  equilibration  in  front  of  seclionA  and  2A  is  not  instantaneous  and  much  of  the  charge  coming 
out  ol  the  feedback  loop,  2A,  stays  at  the  2A  end  and  does  not  divide  properly  Accuracy  is  good  <1% 
error)  but  degrades  with  frequency  and  signal  level.  A  design  improvement  is  described  later  which  im¬ 
proves  the  accuracy  and  speed  of  this  approach  to  levels  similar  to  or  better  than  the  gated  equilibration 
technique. 

The  speed  limitations  of  the  device  are  caused  once  again  by  transfer  inefficiency.  Transfer 
inefficiency  has  been  measured  on  128  stage  CCD  test  devices.  Results  are  shown  in  Figure  9.12  for 
both  the  7.5  fxm  device  and  a  6.0  /jliii  device.  The  6  /ini  device  has  higher  inefficiency  at  low  frequen¬ 
cy  since  the  width  is  also  scaled  down  from  125  /xm  to  100  fxm.  The  speed  improvement  is  slight  and 
indicates  that  further  shrinkage  and/or  buried  channel  processing  are  required  to  attain  the  full  potential 
of  the  charge  domain  concept. 

Due  to  the  small  number  of  stages  required  for  CDF2  the  coefficient  accuracy  holds  up  reasonably 
well  out  to  near  20  MU/.  The  test  bos  has  been  modified  to  run  up  to  35  Mil/  and  the  devices  are  still 
functional  although  the  accuracy  degrades  due  to  transfer  inefficiency. 

Besides  speed  and  accuracy,  the  dynamic  range  is  improved.  Total  dynamic  range  at  I  Mil/  clock 
rate  measures  86  dB  saturation  signal  to  RMS  noise.  In  the  linear  range  of  operation,  80  dB  is  routinely 
achieved. 

9.2  BAND  PASS  FILTER 

following  the  successful  redesign  of  the  low  pass  filter,  a  search  has  been  made  to  find  a  filter  func¬ 
tion  which  can  not  conveniently  be  constructed  with  conventional  technologies  and  pushes  the  limits  of 
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*>.2.1  Killer  Design  <oiixi«lviali»iix 

Both  narrow-  and  wide-hand  multi-pole  bandpass  filters  have  been  investigated  as  potential  test  velu- 
ele  tillers.  The  wide  band  filter  appears  to  be  a  natural  candidate  for  the  charge  domain  approach,  since 
a  single  modified  circular  shift  register  can  provide  six  of  the  desired  poles  simultaneously.  I  he  fre¬ 
quency  response  of  such  a  filter  is  discussed  in  Section  4. 

While  this  approach  yields  an  excellent  filter  with  a  very  efficient  design,  it  does  not  address  the 
generic  case  of  filter  design.  A  narrowband  filler  design,  if  it  can  be  done  in  the  charge  domain,  is  a 
better  challenge  to  test  the  limits  of  CDF  capabilities.  The  task  which  has  been  attempted  is  the  design 
of  a  multi-pole  bandpass  filter  with  a  bandwidth  of  I / 1 00  of  the  clock  frequency  (Fc).  In  addition  the 
passband  itself  should  be  Hat  and  the  skirts  should  be  very  sleep  (down  >  16  dB  at  2/100  l;c 
bandwidth)  to  simulate  the  requirements  of  a  realistic  I'M  communications  system.  The  sidelobes 
should  be  attenuated  more  than  the  40-45  dB  typically  achieved  with  conventional  CCD  split  electrode 
filters  and  the  output  dynamic  range  should  be  kept  as  large  as  possible,  say  60  dB  or  belter.  These 
constraints  turn  out  to  be  quite  difficult  to  meet  with  the  concepts  discussed  thus  far  due  to  some  of  the 
design  considerations  discussed  below. 

fhe  main  advantage  of  using  CDF’s  is  that  all  processing  internal  to  the  filter  is  performed  with 
charge  packet  manipulation.  This  allows  filter  sections  to  be  easily  cascaded  to  build  up  complex  filter 
functions  with  no  overall  loss  in  clocking  speed.  The  drawback  to  this  approach  is  that  because  there 
are  no  gain  elements  in  the  structure,  each  of  the  sections  can  only  attenuate  the  signal.  Thus  if  a 
number  of  stages  are  cascaded  the  attenuation  in  the  passband  (insertion  loss)  may  become  quite  large. 
Since  the  dynamic  range  of  these  devices  is  limited  to  about  80  dB  at  the  input,  this  insertion  loss  can 
severely  degrade  the  filler's  performance. 

I  hese  constraints  have  led  to  a  design  methodology  which  should  be  followed  with  these  devices. 

I  lie  underlying  philosophy  is  that 

when  sections  are  to  be  cascaded,  the  passband  of 
all  sections  should  be  as  similar  as  possible  and 
should  be  designed  for  minimum  insertion  loss  in 
the  pass  hand.  This  will  optimize  the  overall  dynamic 
range. 

I  lie  advantage  of  this  philosophy  is  demonstrated  in  F'-’urc  9.13  which  shows  two  approaches  for  im¬ 
plementing  a  bandpass  filter,  fhe  first  is  a  convention,  cascade  of  single  pole  sections.  In  this  case  the 
A  pole  attenuates  in  the  It  pole’s  passband  and  vice  versa.  Thus  the  overall  passband  has  a  large  inser¬ 
tion  loss  The  second  approach  uses  two  sections  such  as  a  modified  circular  shift  register  cascaded  with 
a  charge  domain  transversal  filler,  where  the  two  sections  have  overlapping  passbands.  The  same 
overall  passband  shape  may  be  achieved  with  much  less  insertion  loss,  thus  significantly  improving  the 
overall  dynamic  range. 

A  second  maior  advantage  of  the  modified  circular  shift  register  can  be  appreciated  in  the  context  of 
tins  design  philosophy,  loi  it  sometimes  permits  more  than  one  of  its  pole  pairs  to  be  used,  thereby  el¬ 
iminating  the  need  to  cascade  two  conflicting  frequency  responses.  I  his  is  accomplished  by  starting 
with  a  circular  shill  register  with  poles  close  to  all  of  the  desired  locations,  and  then  moving  these 
around  In  adding  additional  resonant  channels  so  that  the  pole  locations  simultaneously  move  to  the 
desired  places.  I  he  charge  domain  transversal  filler  is  an  ideal  vehicle  for  providing  many  zeroes  at 
once,  and  the  modified  circular  shift  register  discussed  above  can  be  used  for  implementing  several 
poles  at  once.  I  hese  more  general  types  of  sections  have  enough  flexibility  to  provide  the  desired  fre¬ 
quency  response,  while  greatly  reducing  the  insertion  loss  as  well  as  chip  size  and  power  compared  to  a 
cascade  of  simple  sections. 
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Figure  9.  1 3.  Minimizing  insertion  loss  in  charge  domain  fillers. 

9.2.2  The  Filter  Design  —  CDF3 

The  choice  of  the  narrowband  bandpass  filter,  however,  pushes  beyond  the  limits  ol  this  approach 
and  leads  to  an  extension  of  it,  namely  the  idea  of  using  a  partial  fraction  expansion  as  an  alternative  to 
a  cascade."1 1  In  order  to  make  a  narrowband  filter  with  a  Hat  passband  the  poles  must  be  closely 
spaced.  The  center  frequency  is  chosen  to  he  at  1/4  of  the  clock  frequency  which  leads  to  clustering  of 
the  poles  around  the  imaginary  axis  (or  ()  =  90°  line)  of  the  /-plane.  I  his  choice  yields  the  most  space- 
etlicienl  designs  while  providing  maximum  separation  of  the  alias  band  responses  According  to  the 
equation 

It  =  (2"  1 ) 360/2, V  «  =  1,2,  ■  •  •  ,V  (9.1) 

feedback  loops  with  .V  --  2,  6,  10,  14,  18  .  .  .  stages  of  feedback  result  in  poles  located  at  «  =  90’.  A 
great  deal  of  time  has  been  spent  trying  to  modify  a  circular  shift  register  to  provide  a  close  cluster  of 
poles  around  the  imaginary  axis.  While  this  approach  works  quite  well  for  the  wideband  filler,  the  nar¬ 
rowband  design  calls  for  too  many  stages  in  the  feedback  loop  to  be  practical.  At  this  point  the  fallback 
position  is  again  to  implement  each  pair  of  poles  by  itself  and  cascade  the  sections. 

For  a  bandpass  filler  centered  at  <)  =  90"  and  with  a  bandwidth  of  about  360'7I00  3.6  pole  pairs 

should  be  located  at  approximately  Xl)",  90",  and  91".  The  next  step  is  to  choose  the  amount  of  feed¬ 
back  and  the  number  of  stages  of  feedback  required  lor  each  of  the  sections.  I  he  equation 

(9.31 

implies  that  increasing  the  feedback  portion,  o,  towards  its  maximum  value  of  I  and  increasing  the 
number  of  stages,  V.  will  give  poles  closer  to  the  unit  circle.  However,  increasing  the  feedback  portion 
or  increasing  the  number  of  stages  of  feedback  leads  to  a  larger  chip  si/e  with  higher  power  consump¬ 
tion.  Decreasing  the  (3  on  the  poles  reduces  the  flatness  of  the  passband  and  the  steepness  of  the 
skirls. 

In  order  to  choose  a  "good"  design  a  feedback  limit  of  90%  has  been  chosen  I  he  number  ol 
stages  in  the  feedback  loop  becomes  the  parameter  that  can  be  varied  to  meet  the  performance  goals. 
The  9(1  poles  can  easily  be  implemented  with  ;V  =  2.  so  the  design  of  the  X9  and  91  poles  are  under¬ 
taken  first  A  modified  circular  shilt  register  structure  as  described  in  Section  3  is  used  with  the  majori- 
l\  ol  the  feedback  portion  on  V  \n  2  stages  where  /;  I,  3.  with  a  smallei  portion  on  the 
V  4 in  2  *  I  state,  where  hi  1.2.  and  is  not  necessarib  equal  to  u  I  he  total  feedback 
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Figure  9.  I  7.  Frequency  response  for  cascade  (dashed)  and  parallel  (solid)  connection  of  the  89°  and 
91°  pole  sections. 

must  he  found  to  implement  the  resonances  of  both  modified  circular  shill  registers  in  a  single  struc¬ 
ture.  Although  there  are  several  possibilities  for  doing  this,  a  relatively  simple  solution  has  been  found 
that  proved  satisfactory,  namely  a  parallel  architecture  which  is  equivalent  to  a  partial  fraction  expansion 
of  the  transfer  function.  The  resulting  response,  which  is  compared  to  the  cascade  response  in  Fig¬ 
ure  9.17,  reduces  the  designed-in  insertion  loss  by  13  dB,  which  is  enough  to  make  this  design  practical. 

Two  90'  pole  pair  sections  are  added  in  order  to  flatten  the  passband  and  provide  sideband  attenua¬ 
tion.  The  resulting  parameters  for  these  sections  arc  N  =  2  stages  and  «  =  .875  feedback.  These  sec¬ 
tions  can  be  combined  in  cascade  with  the  89°  and  91°  parallel  combination  since,  as  seen  in  Fig¬ 
ure  9.19,  the  passbands  overlap  in  the  area  of  interest. 

Finally  the  transversal  filter  section  is  added.  The  /.crocs  of  this  filter  arc  placed  on  the  unit  circle 
(see  Figure  9.20)  to  cancel  the  undesired  poles  of  both  of  the  other  sections.  Also,  by  placing  the 
zeroes  symmetrically  around  the  0  =  90°  line,  every  other  coefficient  of  the  transversal  filler  becomes 
zero.  I  he  frequency  responses  of  this  section  (see  Figure  9.21)  also  has  a  passband  which  overlaps  that 
of  the  other  sections  and  can  therefore  be  cascaded  with  little  additional  insertion  loss. 

The  resulting  pole  and  zero  locations  are  shown  in  Figure  9.22  and  are  listed  in  Table  9.23  with  the 
overall  frequency  response  shown  in  Figure  9.24.  Figure  9.25  gives  an  expanded  view  of  the  passband. 
The  original  design  goals  have  been  met  and  the  resulting  filter,  if  it  can  be  fabricated,  promises  im¬ 
pressive  performance.  The  passband  is  flat  to  1  dB  over  1/120  Fc,  the  3  dB  bandwidth  is  1/100  Fc,  the 
response  drops  off  by  18  dB  at  2/100  Fc,  and  by  60  dB  at  5/100  Fc.  The  overall  insertion  loss  is  only 
12  dB  so  that  with  an  input  dynamic  range  of  80  dB,  the  output  stopband  noise  floor  is  -68  dB  with 
respect  to  the  passband. 

It  can  be  seen  from  this  design  exercise  that  the  computer  design  tools  are  invaluable,  not  only  in 
design  optimization,  but  in  the  development  of  the  design  methodology. 


Table  9.23 

POLE  AND  ZERO  LOCATIONS  FOR  THE  BANDPASS  FILTER  DESIGN 


run  cdf4 

Output  file  name:  cdf4 


91  DEG  POLES 

do  you  want  to  include  this?  y 


POLE  LOCATIONS 


# 

REAL 

1 

-0.4352 

3 

0.9684 

6 

0.3621 

7 

-0.9086 

9 

-0.7544 

11 

0.7766 

13 

-0.0328 

15 

0.4141 

17 

-0.7804 

18 

0.0000 

Do  you  want  to  plot  this? 

89  DEG  POLES 

Do  you  want  to  include  this?  y 

POLE  LOCATIONS 

# 

REAL 

1 

0.4245 

3 

-0.9205 

5 

0.0327 

7 

0.9661 

9 

-0.7642 

11 

-0.4212 

13 

0.6255 

15 

-0.2680 

17 

0.7351 

19 

-0.8198 

20 

0.0000 

Do  you  want  to  plot  this? 


IMAG 

RADIUS 

ANGLE 

0.8922 

0.9927 

116.0035 

0.2058 

0.9900 

11.9956 

0.8328 

0.9081 

55.5990 

0.2511 

0.9427 

164.5505 

0.6327 

0.9846 

140.0152 

0.5688 

0.9626 

36.2207 

0.9752 

0.9758 

91.9235 

0.5807 

0.7968 

58.6873 

0.0000 

0.7804 

0.0000 

0.0000 

0.0000 

90.0000 

IMAG 

RADIUS 

ANGLE 

0.8894 

0.9855 

64.4847 

0.2492 

0.9537 

164.8539 

0.9761 

0.9766 

88.0830 

0.2026 

0.9871 

11.8447 

0.6110 

0.9785 

141.3573 

0.8293 

0.9301 

116.9281 

0.4733 

0.7844 

37.1122 

0.7675 

0.8129 

109.2451 

0.5848 

0.9393 

38.5023 

0.0000 

0.8198 

0.0000 

0.0000 

0.0000 

90.0000 

90  DEG  POLES 

Do  you  want  to  include  this?  y 


POLE  LOCATIONS 
# 

1 

2 


REAL  IMAG 

0.0000  0.9354 
0.0000  0.9354 


RADIUS  ANGLE 

0.9354  90.0000 

0.9354  -90.0000 


Do  you  want  to  plot  this? 


90  DEG  POLES 


9-19 


FREQ  (FC) 

0.3222 

0.0333 

0.1847 

0.4571 

0.3889 

0.1006 

0.2553 

0.1630 

0.0000 

0.2500 


FREQ  (FC) 

0.1791 

0.4579 

0.2447 

0.0329 

0.3927 

0.3248 

0.1031 

0.3035 

0.1070 

0.0000 

0.2500 


FREQ  (FC) 

0.2500 

0.2500 


Table  9.23  (Cont’d) 

POLE  AND  ZERO  LOCATIONS  FOR  THE  BANDPASS  FILTER  DESIGN 

Do  you  want  to  include  this?  y 


POLE  LOCATIONS 


# 

REAL 

IMAG 

RADIUS 

ANGLE 

FREQ  (FC) 

1 

0.0000 

0.9354 

0.9354 

90.000 

0.2500 

2 

Do  you  want  to  plot  this? 

0.0000 

-0.9354 

0.9354 

-90.0000 

-0.2500 

ZERO  LOCATIONS 

# 

REAL 

IMAG 

RADIUS 

ANGLE 

FREQ  (FC) 

1 

0.9792 

0.2028 

1 .0000 

11.7000 

0.0325 

3 

0.7785 

0.6277 

1 .0000 

38.8800 

0.1080 

5 

0.4258 

0.9048 

1 .0000 

64.80000 

0.1800 

7 

-0.4258 

0.9048 

1.0000 

115.2000 

0.3200 

9 

-0.7785 

0.6277 

1.0000 

141.1200 

0.3920 

11 

-0.9792 

0.2028 

1.0000 

168.3000 

0.4675 

Do  you  want  to  plot  this? 


POLYNOMIAL  COEFFICIENTS: 


/**() 

C'OEF  - 

1.00000 

/.**  1 

C’OEF  = 

0.00000 

7**2 

COEF  = 

0.98468 

Z**3 

C’OEF  = 

0.00000 

Z**4 

COEF  = 

0.89773 

Z**5 

COEF  « 

0.00000 

Z**6 

COEF  = 

0.97715 

7**1 

COLE  = 

0.00001 

Z**8 

COEF  = 

0.89774 

Z**9 

COEF  = 

0.00000 

Z**I0 

COEF  = 

0.98467 

Z**ll 

COEF  = 

0.00000 

Z**12 

COEF  = 

1 .00000 

TOTAL  RESPONSE 
Do  you  want  to  plot  this? 

9.2.3  Chip  Design  and  Experimental  Results 

The  design  and  layout  of  the  integrated  circuit  chip  which  implemented  the  above  design  is  the  next 
task.  All  of  the  internal  sections  of  the  tiller  have  been  defined.  The  remaining  design  questions  in¬ 
clude  the  signal  input  to  the  chip,  (the  voltagc-lo-charge  conversion)  the  output  chargc-lo-voltage 
conversion,  and  the  connection  of  the  various  sections. 

The  input  section  uses  the  lill-and-spill  technique.  The  fill  operation  is  shown  if  Figure  9.26b  with 
the  excess  charge  spilled  out  in  Figure  9.26c  leaving  a  charge  packet  whose  size  depends  on  the 
difference  in  voltage  between  the  reference  and  input  electrode.  This  charge  packet  is  then  available 
for  clocking  into  the  first  standard  pseudo-one  phase  stage  (Figure  9.26a). 
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Figure  ‘>  26.  Fill  and  spill  charge  input  circuit. 


A  rescl-and-floal  diffusion  is  used  as  an  output  circuit  with  a  two  stage  source  follower  to  drive  olT- 
cliip  electronics  or  scope  probe  capacitance.  The  schematic  is  shown  in  Figure  9.27. 

I  he  usual  technique  of  connecting  isolated  sections  of  a  CCD  chip,  namely  using  a  diffusion  at  the 
inputs  and  outputs  ol  the  subsections  connected  together  with  a  metal  line,  gives  errors  in  the  real  pole 
section  ol  CDF2.  These  errors  are  caused  by  a  lack  of  lateral  equilibration  of  the  charge  packet  prior  to 
the  harrier  splitting  operation.  For  a  pseudo  one  phase  CCD  clocking  scheme  there  arc  two  DC  gates,  a 
iranslei  gale  and  a  storage  gate.  Father  or  both  of  these  may  be  sliced  into  two  halves  in  the  direction 
perpendicular  to  charge  How  in  the  channel,  and  a  dilFused  n-lype  region  may  be  placed  in  between  the 
tun  halves  in  order  to  enhance  the  lateral  equilibration  speed.  The  dillusion  may  also  have  a  metal 
overlaying  n  to  lurther  enhance  litis  speed  Flits  is  shown  in  Figure  9.28  where  both  the  DC  transfer 
and  DC  storage  gates  are  split.  A  metali/ed  dillusion  may  also  be  placed  between  the  transfer  and 
storage  gates,  without  significantly  degrading  the  transfer  speed.  The  metallized  dillusion  can  actually 
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be  iwo  di fl used  areas  connected  by  the  metal,  as  is  done  in  turnarounds,  etc.  In  previous  designs  such 
as  CDF2  a  single  metallized  diffusion  is  used  as  the  interconnection  of  lion-adjacent  channels  as  well  as 
a  mechanism  lor  high  speed  lateral  equilibration.  This  causes  the  _*edback  coefficient  to  be  in  error. 


The  new  structure  which  splits  both  the  transfer  and  storage  electrodes  separates  the  functions  of 
connection  and  lateral  equilibration  1  he  connection  function  is  performed  by  the  split  transfer  elec¬ 
trode.  T  he  charge  then  spills  into  the  split  storage  electrode  so  that  even  if  all  of  the  charge  is 
transferred  to  one  end  of  the  storage  electrode,  once  there  it  will  equilibrate  along  its  whole  length  be¬ 
fore  being  transferred  on  to  the  next  stage  where  barrier  splitting  occurs  This  type  of  interconnect 
avoids  the  worries  of  balancing  capacitance  and  minimizing  resistance.  This  technique  is  used  for  con¬ 
necting  all  the  sections  ol  the  bandpass  liltei  In  addition  this  slim  line  has  been  lesled  as  a  new  imple 
nielli,  illoll  ol  (tie  <  I  >  I  .’  leal  pole  Willi  hall  lei  .plllllll)'  and  silii'le  slape  Ireilb.n  t  l-.ei  I  igtlle  '*  "•>  I  In 
n  i  in  ai  \  ol  the  leal  puli  i  m  I  Id  lent  i  .  iinpi  oinl  1 1  mil  ei  i  m  In  .1  II  1  mu  and  -.1  low  s  tu  1  sir  lillli  ml 
Mill  main  III  willi  I  n.  is  le  \  el  ni  i  li  n  k  speed  (  up  to  I  lie  Inn  Ms  nil  pi  >sed  In  11  ansi  ei  me  I  hi  lein  \  > 
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Figure  *>. 2S.  CCD  structure  Cor  enhancing  lateral  equilibration. 


The  chip  (1)1-3  has  been  layerl  out  according  to  the  design  considerations  discussed  above.  The 
chip  (Figure  9.30)  has  been  processed  with  the  (>  pun  design  rule  n-channel  process  as  used  for  CDF2. 
An  additional  test  chip  is  included  to  allow  testing  of  the  individual  filter  sections.  The  test  chip  demon¬ 
strates  proper  operation  of  each  of  the  sections  individually,  verifying  operation  of  the  modified  circular 
shill  register  as  well  as  the  new  method  of  interconnecting  filler  sections.  The  impulse  response  mea¬ 
surements  again  show  accuracy  to  about  0.2%  anti  the  frequency  response  as  measured  on  a  spectrum 
analwcr  agrees  with  the  simulation  to  belter  than  0.5  dll  from  DC  up  to  the  Nyquist  frequency. 

I  he  overall  filter,  however,  does  not  operate,  flic  problem  is  due  to  inadequate  handling  of  the 
bias  charge.  While  the  early  stages  may  attenuate  the  signal  charge  level,  the  bias  charge  is  not  at¬ 
tenuated.  l  ater  stages  must  be  built  to  handle  the  entire  50%  of  lull  well  (FW)  bias  charge  of  the  input 
stage.  I  he  overall  charge  level  at  the  output  is 

-  .5  FW  ±  (.5  x  SAIT 

where  SAI  is  the  signal  attenuation  factor  in  the  passband. 

The  circuit  has  been  redesigned  as  CDF4  with  several  improvements.  The  chip  photo  is  shown  in 
Figure  9.31.  The  design  has  been  converted  to  4  pun  design  rules  which  gale  a  transfer  gate  length  of 
4  pun  and  storage  gale  length  of  9  pun  versus  6  and  10  p.m  for  CDF3.  The  stages  which  give  no  inser¬ 
tion  loss  (the  transversal  filler  and  90"  pole  sections)  are  moved  to  the  front  end  in  order  to  improve 
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Figure  9.31.  CDF4  bandpass  filter  photomicrograph. 

their  accuracy.  With  this  placement  scheme,  the  device  handles  a  greater  ratio  of  signal  charge  to  bias 
charge,  and  bias  charge  variations  are  not  as  critical.  The  input  is  made  directly  into  the  transversal 
filler  section.  The  channel  width  is  320  (im,  which  is  the  minimum  channel  width  which  gives  good  ac¬ 
curacy  in  the  splitting  operation.  The  next  section  is  a  90°  pole  pair  section.  The  overall  channel  width 
is  1280/xw  with  a  feedback  coefficient  of  0.875.  The  steady  state  bias  handling  channel  has  a  width  of 
(0.125)  1280  fxiu  =  160  fxm.  Thus  interposed  between  the  transversal  filter  section  and  the  90°  pole 
section  is  a  section  which  divides  the  packet  into  two  halves  and  throws  one  half  away.  Following  the 
first  90°  pole  section  is  a  second  90°  pole  section.  Its  output  is  split  in  halves  with  one  half  directed  to 
the  89°  pole  section  and  the  other  half  connected  to  the  91°  pole  section.  The  channel  widths  of  these 
two  sections  are  800  fim  with  a  total  feedback  coefficient  of  0.9.  The  bias  handling  channel  is 
0.1(800)  =  80  ftm  wide.  These  two  sections  have  cross-connected  inputs  to  provide  the  subtraction 
operation.  The  outputs  are  summed  together  and  fed  to  the  output  amplifier,  a  reset-and-float  structure 
with  a  two  stage  source  follower. 

These  redesigned  chips  have  been  fabricated  in  General  Electric’s  Corporate  Research  and 
Development’s  semiconductor  processing  facility.  The  frequency  response  of  the  CDF4  chips  as  mea¬ 
sured  on  an  HP8553A  spectrum  analyzer  is  shown  in  Figure  9.32.  Figure  9.33  compares  this  response 
to  the  simulation.  The  agreement  is  excellent,  with  errors  less  tljan  1.0  dB  outside  the  noise  limited  re¬ 
gime.  The  limit  on  stopband  attenuation  of  65  dB  is  set  by  device  noise  rather  than  by  coefficient  accu¬ 
racy.  Measured  insertion  loss  is  12.0  dB  versus  the  12.1  dB  predicted.  An  expanded  view  of  the 
passband  (f  igure  9.34)  shows  good  agreement  between  experimental  and  theoretical.  The  experimental 
results  show  some  extra  peaking  in  the  89°  pole  compared  to  the  91°  pole.  This  is  caused  by  a  uneven 
division  of  the  charge  which  is  sent  to  the  89°  and  91°  pole  sections.  Figure  9.35  shows  the  effect  of  a 
0.5%  error  in  the  division  into  two  halves.  A  redesigned  chip  should  make  this  split  separately  from  the 
1/8  —  7/8  split  used  for  the  90°  pole  in  order  to  make  the  50%  split  completely  symmetric. 
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Figure  9.32.  CDF4  experimental  frequency  response. 

The  demonstration  of  this  filter  dramatically  shows  the  unique  capabilities  of  charge  domain  devices 
for  monolithic  signal  processing.  In  addition  to  demonstrating  exceptional  frequency  response  charac¬ 
teristics,  this  filter  verified  our  simulation  tools  as  well  as  our  capability  to  determine,  in  advance,  the 
filter  characteristics  which  can  be  achieved  with  the  charge  domain  technique. 
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Figure  9.33.  Comparison  of  theoretical  (solid)  and  experimental  (dotted)  frequency  responses  for 
CDF4. 


Figure  9.34.  Comparison  of  theoretical  (solid)  and  experimental  (dotted)  frequency  responses  for 
CDF4  —  expanded  view. 


9-28 


Fraquenoy  (Fo) 


i 

► 

i 

j 

» 

► 

r 

i 


I  igurc  9.35.  The  e  tie  el  of  unequal  charge  packets  transferred  to  the  89°  ami  *>|  pole  sections  on  the 
passhanil  of  (Dl'4  for  0.0%  error  (solid  curve)  and  0.5""  error  (dashed  curve). 
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10.  PROGRAMMABLE  CDFS 

The  filters  described  in  the  previous  section  show  good  performance,  but  are  limited  in  their  flexibil¬ 
ity  or  adaptability  since  their  coefficients  and  configurations  are  fixed  in  their  design.  It  is  desirable  in 
many  cases  to  electrically  program  the  filter  coefficients  or  change  the  filter  response  after  it  has  been 
fabricated.  A  number  of  techniques  have  been  investigated  in  order  to  achieve  this  goal.  Programma¬ 
ble  test  devices  have  been  fabricated  along  with  the  CDF3  and  CDF4  filters  (see  Figure  10.1a).  The 
evaluation  of  these  devices  proves  that  the  concepts  work  as  expected.  It  now  appears  practical  to  make 
a  filter  which  combines  the  benefits  of  charge  domain  signal  processing  with  the  flexibility  of  electrically 
programming  the  coefficients.  This  accomplishment  coupled  with  recent  developments  in  microcom¬ 
puters  and  microcontrollers  opens  up  a  new  realm  of  possible  applications  for  charge  domain  signal  pro¬ 
cessing  chips. 

10.1  ALTERNATE  PROGRAMMABILITY  TECHNIQUES 

There  are  several  degrees  of  programmability  which  can  be  incorporated  in  a  CDF  design,  and  in 
general  the  more  flexible  approaches  are  more  costly  in  terms  of  chip  size,  speed,  or  complexity.  The 
lowest  degree  of  programmability  is  inherent  in  any  clocked  CDF  design.  The  absolute  frequency 
response  can  be  altered  just  by  changing  the  clocking  frequency.  A  bandpass  filter  can  shift  its  center 
frequency  in  proportion  to  the  clock  frequency.  The  bandwidth,  however,  also  scales  with  the  clock 
frequency.  For  some  communications  systems,  where  the  tuning  range  is  small,  this  degree  of  flexibility 
may  be  adequate. 


Figure  10.1a.  Charge  domain  chips  including  CDF3  filter  and  PCDF  programmable  test  chips. 
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The  next  step  up  in  flexibility  is  to  have  a  fixed  band  shape  but  have  independently  variable  center 
frequency  and  bandwidth.  An  approach  which  implements  this  function  is  the  N-path  filter  shown 
schematically  in  Figure  10.2.  The  commutator  at  the  input  demodulates  the  desired  incoming  signal 
frequency  down  to  baseband  samples  and  the  commutation  speed  determines  the  center  frequency. 
The  baseband  samples  pass  through  a  bank  of  parallel  low  pass  charge  domain  filters.  The  clocking  fre¬ 
quency  of  these  filters  is  used  to  program  the  bandwidth.  The  output  of  these  filters  can  be  synchro¬ 
nously  modulated  up  to  the  original  frequency  if  desired.  Alternatively  the  baseband  signals  may  be 
used  as  outputs.  When  the  number  of  parallel  filters  is  4,  the  baseband  signals  can  be  combined  to  give 
the  in  phase  (I)  and  quadrature  (Q)  outputs  directly.  This  technique  has  been  analyzed  to  some  extent, 
but  most  of  the  efforts  in  this  area  have  been  directed  at  the  next  tier  of  programmability. 

The  next  level  up  in  the  programmability  chain  is  the  capability  of  altering  the  coefficients  electrical¬ 
ly,  but  not  changing  the  architecture.  The  first  task  is  to  define  programmable  charge  splitter  struc¬ 
tures.  Three  basic  structures  have  been  invented  and  evaluated  experimentally.  The  first  is  an  in-place 
sequential  charge  splitter.  This  approach  uses  a  gated  equilibration  splitting  scheme  similar  to  that  used 
successfully  in  CDF2.  A  digital  shift  register  supplies  the  coefficients,  most  significant  bit  (MSB)  first. 
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Figure  10.1b.  PCDF1A  =  in-phase  sequential  splitter 

PC'DFIB  —  four  stage  pipe  organ  filter  using  sequential  splitters. 
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A  schematic  diagram  of  this  technique  is  shown  in  Figure  10.3.  The  charge  packet  is  introduced  into 
potential  well  under  the  A  gate.  Meanwhile  the  B  well  is  empty.  At  this  time  the  equilibration  gate  is 
opened  and  closed  leaving  half  of  the  charge  packet  under  A  and  half  under  B  (note  that  a  splitting  ra¬ 
tio  of  0.5  simplifies  the  balancing  of  the  charge  splitting  mechanism).  At  this  point  either  the  C  or  D 
gate  is  turned  on  depending  on  the  (MSB),  and  the  1/2  charge  packet  is  transferred  to  the  E  or  F  accu¬ 
mulator  well.  Now  the  A  well  contains  1/2  the  input  packet  and  the  B  well  is  empty.  The  equilibration 
gate  again  turns  on  and  off  leaving  1/4  of  the  original  packet  under  A  and  1/4  under  B.  The  second  bit 
of  the  digital  coefficient  channels  the  1/4  under  B  to  the  E  or  F  accumulator.  The  splitting  and  chan¬ 
neling  continues  for  1/8,  1/16,  etc.,  for  /V  bits  where  N  is  as  many  bits  as  is  desired.  After  all  the  sig¬ 
nal  (possibly  including  the  leftover  1/2^  portion)  has  been  channeled  the  output  can  be  read  out  with 
an  overlying  electrode  (non-destructive  readout)  and  the  entire  charge  packet  returned  to  the  A/B 
reservoir,  or  with  a  diode  (destructive  readout)  and  a  new  input  sample  can  be  introduced.  This  tech¬ 
nique  has  the  advantages  of  being  compact  in  size  and  universal  in  structure  with  respect  to  the  number 
of  bits  that  may  be  used.  Its  disadvantage  is  that  the  cycle  time  for  a  full  multiplication  is  equal  to  the 
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Figure  10. Ic.  PCDF1C  =  four  stage  split  electrode  filter  using  sequential  splitters 
PCDF1D  =  flash  splitter 
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clock  period  of  the  splitter  multiplied  by  the  number  of  bits.  This  disadvantage  is,  however,  ameliorat¬ 
ed  by  the  high  inherent  speed  of  the  splitting  mechanism. 


One  of  the  test  devices,  PCDF1A  (Figure  10.1b)  has  been  used  to  evaluate  this  technique.  For  this 
and  all  other  PCDF  test  chips  the  digital  shift  registers  are  off-chip  although  a  production  circuit  would 
include  on-chip  shift  registers.  The  device  is  operated  as  a  multiplying  digital-to-analog  converter 
(MDAC)  where  the  analog  input  is  a  DC  value  and  the  coefficients  are  decremented  from  127  to  0  (for 
a  7  bit  conversion).  The  output  shown  in  Figure  10.4  demonstrates  excellent  linearity  over  the  entire 
range,  with  the  greatest  error  at  the  64-63  coefficient  transition.  Seven  to  eight  bits  of  linear  operation 
have  been  demonstrated  with  no  degradation  at  splitting  speeds  up  to  15  MHz.  The  number  of  bits  is 
limited  by  the  accurate  transfer  of  charge  from  the  B  well  to  the  E  or  F  wells  in  a  short  period  of  time. 
Structural  improvements  can  yield  9-10  bits  of  capability. 


The  second  technique  overcomes  the  speed  limitations  of  the  sequential  splitter  at  the  expense  of  an 
increase  in  chip  size.  This  is  the  pipelined  charge  splitter  technique  diagramed  in  Figure  10.5.  In  this 
case  splitting  takes  place  unidirectionally  in  several  stages  with  each  stage  providing  one  binary  bit  of  in- 
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creased  resolution.  The  A  and  A  accumulators  are  clocked  in  synchronism  with  the  splitters  such  that 
after  8  clock  cycles  the  product  result  appears  at  the  output.  Due  to  pipelining  the  next  result  appears 
one  cycle  later.  Therefore  the  throughput  rate  of  this  technique  is  equal  to  the  clock  rate. 

The  third  approach  is  a  flash  splitter.  In  this  approach,  shown  in  Figure  10.6,  the  channel  is  simulta¬ 
neously  split  into  portions  equal  to  1/2,  1/4,  1/8,  .  .  .  1/2  N  and  another  1/2  N  times  the  input  chan¬ 
nel.  This  approach  gives  a  full  n-bit  multiplication  in  one  clock  cycle.  It  requires  a  very  wide  structure 
in  order  to  accurately  divide  the  charge. 

Both  of  these  concepts  have  been  tested  out  experimentally.  The  single  bit  per  stage  device, 
PCDF1E  (Figure  10. Id)  shows  operation  as  expected,  again  with  7-8  bits  of  accuracy  in  steady  state 
tests.  The  device  is  operated  as  an  MDAC  with  a  DC  analog  input  and  varying  digital  coefficients.  Un¬ 
like  the  sequential  splitter,  the  accuracy  is  much  poorer  than  7-8  bits  if  the  digital  coefficient  changes 
from  one  cycle  to  the  next.  This  error  is  caused  by  the  lack  of  introduction  of  a  fat  zero  charge  into  the 
running  accumulators.  The  high  transfer  inefficiency  results  in  a  large  error  for  the  first  output  value 
after  a  coefficient  changes.  Figure  10.7  shows  the  results  for  a  constant  DC  input  and  eight  successive 
outputs  (both  plus  and  minus  outputs)  for  a  coefficient  of  0,  eight  outputs  for  a  coefficient  of  1 ,  and  so 
on  up  to  a  coefficient  of  15.  Note  the  error  in  the  first  output^ for  each  group.  This  problem  is  due  to 
high  charge  transfer  inefficiency  in  the  initially  empty  A  and  A  accumulators.  This  problem  can  easily 
be  alleviated  in  future  designs  by  adding  a  fat  zero  charge  into  both  accumulators.  The  four  bit  flash 
design  PCDF1D  (Figure  10.1c)  has  also  been  evaluated.  It  exhibits  6-7  bits  of  linearity  in  the  MDAC 
mode  (see  Figure  10.8).  The  accuracy  appears  to  be  limited  by  cross  coupling  of  the  address  lines  into 
the  output.  Careful  shielding  and  on-chip  digital  shift  registers  should  alleviate  this  problem.  No 
transfer  inefficiency  effects  are  apparent  with  this  design. 
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Figure  10.6.  Flash  charge  splitter  schematic. 
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Figure  10.8.  Flash  charge  splitter  operated  as  a  four  bit  MDAC. 


The  optimum  choice  for  a  practical  chip  may  be  a  hybrid  combination  of  the  pipeline  and  flash  ap¬ 
proaches.  A  compact,  high  speed  9  bit  multiplier  cell  may  include  three  pipelined  stages  with  three  bits 
of  flash  conversion  per  stage.  This  approach  appears  to  be  quite  practical  if  careful  shielding  and  proper 
bias  charges  are  utilized. 

The  next  concern  in  a  circuit  implementation  is  configuration,  or  how  should  these  programmable 
coefficient  blocks  be  hooked  up.  Two  representative  approaches  have  been  tested.  Figure  10.9  is  a 
schematic  of  a  4-channel  pipe  organ  type  connection12’  such  as  is  required  for  a  charge  domain  pro¬ 
grammable  transversal  filter  or  charge  domain  recursive  filter  section.  The  test  circuit  PCDF1B 
(Figure  10.1b)  is  such  a  structure  which  uses  sequential  splitters  to  implement  the  programmable 
coefficients.  The  output  of  the  four  channels  is  summed  and  sensed  with  a  resel-and-float  diffused 
diode.  Each  cell  operated  as  described  previously  with  7-8  bits  of  linearity  with  the  output  equal  to  the 
sum  of  the  outputs  from  each  section.  The  accuracy  is  good  although  again  better  shielding  should  be 
incorporated. 
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figure  10.9.  Pipe  organ  type  programmable  transversal  lilter. 

The  alternate  circuit  configuration,  shown  in  f  igure  10.10  is  similar  to  a  conventional  CCI)  transver¬ 
sal  tiller.  In  this  ease,  a  sequential  splitter  is  used  to  give  a  programmable  tap  weight.  As  in  the  con¬ 
ventional  split  electrode  filler,  the  charge  is  sensed  non-destruclively  on  overlying  electrodes.  Once  the 
charge  packet  is  sensed,  it  may  be  transferred  on  to  the  next  stage.  Alternatively,  the  charge  packet 
mat  be  retained  at  one  site  and  the  digital  coefficients  may  be  shilled  from  one  stage  to  the  next,  the 
so-called  moving  reference  structure.  The  first  approach  is  included  as  a  lest  chip.  PC  Dl  1(  . 
(figure  10. 1  c >  and  operates  as  expected,  the  output  being  the  sum  of  the  outputs  of  the  individual 
cells.  I  lie  second  approach  has  some  inherent  advantages  over  this.  I  list  of  all.  charge  transfer 
inefficiency  effects  will  not  degrade  performance  since  the  charge  packet  is  retained  in  a  single  cell. 
Secondly,  the  loading  and  shifting  of  the  digital  coefficient  is  easier.  Third,  and  most  important,  is  the 
ability  to  easily  connect  up  many  sections  and  even  many  chips  in  series  to  obtain  very  long  correlation 
products  or  huge  time-bandwidth  products.  One  such  design  is  discussed  in  the  next  section. 
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I  igure  10.10.  Splii  electrode  type  programmable  transversal  filler 


II.  APPLICATIONS 


This  contract  lias  addressed  the  understanding  and  development  of  tools  lor  designing  (  I >1  s.  |  he 
next  task  is  to  find  applications  which  could  utilize  charge  domain  devices  to  icduce  system  size,  powci. 
and  cost  and/or  improve  performance.  Several  features  of  ('l)l  \  which  are  important  loi  some  appli¬ 
cations  are  discussed,  and  then  several  representative  applications  ate  detailed. 

II. I  FURTHER  CAPABILITIES  OF  C  DF'S 

In  the  drive  for  higher  speed  monolithic  signal  processing,  several  features  of  charge  domain  devices 
become  particularly  attractive.  The  first  advantage  which  can  be  utilized  is  the  higher  speed  and  lowei 
transfer  inefficiency  of  buried  channel  device  fabrication.  As  discussed  earlier,  linearity  is  maintained 
with  buried  channel  CDF's  and  their  clocking  rate  is  only  limited  by  transfer  efficiency  effects. 

Another  feature  of  CDF's  in  particular  and  sampled  data  systems  in  general  is  the  use  of  alias  band 
processing,  where  the  incoming  signal  frequency  is  above  the  device  clocking  frequency.  For  instance 
the  bandpass  filter  CDF4  has  its  center  frequency  </,,)  at  1/4  of  the  clock  frequency  f/L  However 
frequencies  of  3/4  /,,  5/4  /, ,  7/4  /,  or 


./  =  «/',  3- 


n  -  0,1,... 


will  also  be  passed.  In  addition  the  sampling  operation  can  be  used  to  down-convcn  the  signal  to  the 
n  =  0  band.  Short  prefillers  can  be  used  to  suppress  the  undesired  passbands 

11.2  APPLICATIONS  EXAMPLES 

One  application  which  can  use  this  technique  is  shown  in  Figure  ILL  The  primary  filter  which  does 
the  band  shaping  is  the  75  stage  filter  which  is  clocked  til  20  MHz.  It  has  its  passbands  at  /  =  5,  15, 
25.  .  .  .  MHz.  The  CDF4  chip  can  be  used  as  this  filler.  Short  prelilters  can  be  added  on  the  same  chip 
to  suppress  those  passbands  at  5,  15,  .  .  .  65,  K5,  95,  .  .  .  i.c.,  just  leaving  the  75  MHz  passband.  The 
prelilters  only  require  a  few  charge  domain  transversal  filter  stages  and  therefore  require  little  power 
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MAIN  FILTER:  75  STAGES  -  CLOCKED  AT  20  MHZ 


SECOND  PREFILTER:  5  STAGES  -  CLOCKED  AT  80  MHZ 


THIRD  PREFILTER  5  STAGES  -  CLOCKED  AT  40  MHZ 


OVERALL  RESPONSE  -  2  MHZ  BANDWIDTH  AT  75  MHZ 


I  igurc  I  1.1  75  Mllz  il  filtering  utilizing  alias  bands  and  small  prelilters. 


iimI  in.  lelalively  insensitive  In  lianslei  mclliciency  |{y  using  alms  hands  and  prelillcrs  a  high  pcilor- 
iiiiinu'  luw  power  charge  domain  fillet  can  he  marie  to  operate  at  very  lugh  speeds. 

Another  application  of  this  technology  is  to  improve  performance  of  conventional  CCD  split  elec¬ 
trode  filters.  The  required  zeroes  of  a  transversal  filler  can  be  divided  between  the  charge  domain  sec¬ 
tion  and  the  split  electrode  section.  Such  a  structure  is  diagrammed  in  Figure  11.2.  Oftentimes  the 
Idler  function  can  he  factored  such  that  the  CDF  section  requires  only  positive  coefficients,  further  sim¬ 
plifying  the  operation.  One  practical  application  of  this  technique  is  to  use  the  higher  speed  CDI  as  an 
anti-alias  preli Iter  (as  in  the  last  example)  with  the  lower  speed  split  electrode  filter  or  programmable 
CDF  transversal  filter  providing  the  band  shaping  function. 

Another  application  which  uses  the  alias  band  processing  as  well  as  the  demodulation  capabilities  of 
the  sampling  operation  is  an  FM  communications  system  chip.  This  application  has  been  analyzed  re¬ 
cently  and  a  white  paper  tiescribing  the  system  was  submitted  in  March  1983.  This  paper  is  included 
here  as  Appendix  4 

Finally  a  programmable  charge  domain  device  has  been  analyzed  using  the  moving  reference  archi¬ 
tecture  with  sequential  splitters.  One  hundred  stages  can  be  incorporated  on  a  chip.  This  chip  can 
function  as  a  programmable  digital-analog  correlator,  a  programmable  or  adaptive  finite  impulse 
response  tiller,  or  a  pulse  compression  filter.  Several  of  these  chips  can  be  cascaded  for  longer  correla¬ 
tion  products.  This  chip  was  described  in  a  white  paper  report  submitted  January  1983  and  is  included 
here  as  Appendix  5. 

f his  list  of  applications  is  by  no  means  all  inclusive;  it  is  only  a  suggestion  of  what  may  now  be 
practical  with  the  demonstration  of  the  capabilities  of  charge  domain.  Other  application  areas  for  which 
this  technology  may  be  attractive  include  IF  filters,  21)  programmable  correlators  and  image  processors, 
steerable  radar  or  sonar  beam  formers,  spread  spectrum  communication  coders  and  decoders,  video 
filters.  FIT's  matrix  multipliers,  D/A  converters,  A/D  converters  or  any  function  requiring  high  data 
rale  sampled  data  signal  processing. 
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Figure  11.2.  CDI  /(’CD  hybrid  device. 


12.  SUMMARY 


The  contract  effort  to  develop  charge  domain  filters  lor  enhanced  monolithic  signal  processing  was 
quite  successful.  The  first  goal  of  investigation  ol  the  capabilities  and  limitations  of  these  devices  was 
addressed  through  studies  of  the  charge  splitting  operation,  studies  of  various  specific  charge  manipula¬ 
tion  structures,  studies  of  the  use  of  poles  as  well  as  zeroes  for  implementing  transfer  functions,  and 
studies  of  techniques  for  placing  poles  and  zeroes  in  the  /-plane  efficiently  under  the  constraints  of  the 
charge  domain  technology.  In  addition  testing  techniques  and  simulation  tools  were  developed  10  aid  m 
the  design  and  understanding  of  operation  of  charge  domain  filters. 

The  second  goal  of  the  program  was  the  actual  design,  lav  out,  ami  fabrication  of  a  multi-pole  liltct 
which  would  demonstrate  the  unique  capabilities  of  charge  domain  technology.  A  narrowband  bandpass 
filter  was  investigated  as  the  lest  vehicle,  since  it  pushed  beyond  the  limits  of  what  was  through  pravti 
cal  at  that  point.  Several  key  discoveries  permitted  this  filter  to  be  designed  with  impress, vc 
specifications.  The  device  was  fabricated  and  the  demonstrated  performance  matched  the  theotetie.il 
predictions  to  a  high  level  of  accuracy.  A  general  design  methodology  for  ('l)T’s  was  developed  and 
expanded  which  along  with  the  simulation  tools  permits  application  to  be  investigated  quickly  and  prac¬ 
tical  chips  to  be  designed  easily. 

The  final  goal  of  the  contract  was  the  investigation  of  techniques  which  would  enable  a  C ' I ) I  to  be 
programmed  electrically,  presumably  by  a  microcomputer  in  a  systems  application.  A  number  of  levels 
of  programming  flexibility  were  investigated  and  some  concepts  were  designed  in  silicon  and  tested  out. 
These  programmable  coefficient  test  blocks  and  four  stage  filters  successfully  demonstrated  the  ability  to 
provide  7-8  bits  of  digital  coefficient  accuracy  in  the  multiplication  operation.  Large  multi-stage  chips 
with  8-10  bits  of  accuracy  appeared  to  be  practical. 

f  inally  a  number  of  applications  areas  and  several  specific  applications  were  investigated  as  potential 
follow  on  programs  to  utilize  these  capabilities.  These  ranged  from  I  'M  fillers  and  demodulators  to  low 
power  video  filters  to  programmable  II)  and  21)  correlators/lillers.  T  he  potential  for  incorporating  high 
performance  high  speed  systems  on  a  chip  utilizing  charge  domain  technology  is  enormous. 
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Appendix  1 

TIME  DOMAIN  TYPE  TEST  PROGRAM 

10  PRINT:PRINT 
15  DEFINl  A-N:DEPDBL  S-T 
18  GOTO  1000 

20  INPUT  2WHICH  PROGRAM  81  =CDF2,2  =  7)”*N 
30  ON  N  GOTO  1000,1000 
40  GOTO  20 
1000  REM 

1010  REM  CDF1  10/26/81 
1020  REM 

1030  REM  THIS  PROGRAM  CALCULATES  AND  PRINTS 

1040  REM  THE  IMPULSE  RESPONSE  OF  CDF2 

1050  REM 

1060  NS=  128 

1070  DIM  X(NS),Y (NS) 

1110  INPUT  “#READS,  #  D1SPLAYED”;NR,NP 
1140  SCALE  =  5/4096 
1150  PRINT 

1160  PRINT  “  DATA  ”,  “OUTPUT”, “NOISE”,  “OUTPUT,"  NOISE  ” 

1170  PRINT  “POINT  #“,”  VOLTS”, “VOLTS”, “MAXOUT=  l”,“MAXOUT=  1” 

1170  PRINT  “ _ ”,“ _ ”,“ _ ”,“ _ ”,“ _ ” 

1180  MM  =  114000 

1183  REM  GET  DATA  SUMS 

1190  FOR  1=1  TO  NS 

1190  M  =  H9000  +  I-1 

1193  X A=O:X3  =  0 

1200  FOR  J  =  1  TO  NR 

1220  A=16*PEE<(M)  +  PEEK(M  +  MM) 

1222  XA  =  XA  + A 

1224  XB  =  BB  +  Z*A 

1225  4  =  M  +  NS 
1230  NEXT 

1240  X(I)  =  SCALE*X A/NR 

1241  SD  =  (XB-XA*XA/NR)/(NR-1) 

1242  Y(I)=SCALE*SQR(ABS(SD)) 

1243  OUT  4,1 
1250  NEXT 

1255  REM  SUBTRACT  ZERO  AND  FIND  MAX 

1260  SX=0:A=0 

1270  FOR  I  =  3*NX/4  +  I  TO  NS 

1280  SX  =SX*X(I) 

1290  NEXT 

1300  SX=SX/(NS/4) 

1310  FOR  1  =  1  TO  NO 
1320  Z(I)  =  X(!)-SX 

1340  IF  ABS(X(1))> YM  THEN  YM=ABS(X(I)> 

1350  NEXT 

1360  IF  YZ  =  0  THEN  YM  =  100 
1370  FOR  1  =  1  TO  NP 

1380  PRINT  1,<(I),Y(I),X(I)/YM,Y(I)/YM 
1390  NEXT 
1390  PRINT 

1390  AS  =  “Y”:BS  =  “Y” 
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Appendix  2 

DIGITAL  SIGNAL  PROCESSNG  ROUTINES 

C  DSP  4/11/83 
C 

C  THIS  PROGRAM  CALCULATES  AND  PRINTS  AND/OR  PLOTS 
C  THE  IMPULSE  RESPONSE, FREQUENCY  RESPONSE, AND 
C  POLE-ZERO  LOCATIONS  FOR  A  CDF 
C  GIVEN  THE  POLE-ZERO  LOCATIONS  OR 
C  THE  TRANSFER  FUNCTION  H(Z) 

C 

C  SUBROUTINE  LOCATIONS: 

C  IMSL  SBR-ZPOLR(COEF,NDEG, ROOTS, 1ER) 

C  1500  SBR-PZIMP(ZERO,NZERO) 

C  1800  SBR-FREQPZ(ZFREQ,NFREQ,ZERO,NZERO) 

C  2000  SBR-PZ(POLE,NPOLE,ZERO,NZERO,PCOEF,N,H) 

C  2500  SBR-COEF(PCOEF,NPOLE,ZCOEF,NZERO,N,H) 

C  3000  SBR-M  AG(X,XMAG,XMAGDB,PH1,N) 

C  3500  SBR-SAMPLEIXM AG,XM  AGDB,PHI,N) 

C  4000  SBR-FFT(X.N) 

C  4500  SBR-PZMINMAX(POLE,NPOLE,ZERO,NZERO,PCOEF,N,Y) 

C  5000  SBR-PRNTIMP(XTIME,NT) 

C  5300  SBR-PRNTFREQ(XM  AG,XM  AGDB,PHI,N) 

C  5600  SBR-PRNTPZ(POLE,NPOLE,ZERO,NZERO) 

C  6000  SBR-PLTIMP(XTIME,NT,KLAB, SCALE) 

C  6300  SBR-PLTMAGIXM  AGDB,N,KLAB,SCALE) 

C  7000  SBR-PLTPHASE(PHI,N,KLAB, SCALE) 

C  7500  SBR-PLT(ARRAY,N,KLAB,SCALE,XMIN,XMAX,XINC, 

C  YMIN,YMAX,YNC,XPOINT,YPOINT,NPOINT,KCH  AR) 

C  NOTE:  N  =  1  FOR  IMP;  KLAB  =  NEG  FOR  PHASE 
C  NOTE:  SCALE=  1.  FOR  8X10  PLOT;SCALE  =  -l .  FOR  16X10  PLOT 
C  KPEN  PICKS  THE  PEN  COLOR 

C  KLAB-THOU:  0  =  ASKS  FOR  TITLE  /  1=  AUTOMATIC  PLOT 
C  KLAB-HUND:  0  =  NO  AXES  /  KPEN  =  AXES  AND  LABELS 
C  KLAB-TENS:  0  =  NO  CURVE  /  KPEN  =  PLOT  CURVE 

C  KLAB-ONES:  0  =  NO  POINTS  /  KPEN  =  PLOT  DATA  POINTS 

C  9000  SBR-PLTPZ(POLE,NPOLE,ZERO,NZERO.KLAB,SCALE) 

C 

****************************************** 

C 

SUBROUTINE  PZIMP(ZERO,NZERO) 

C  THIS  ROUTINE  CONVERTS  THE  ZERO  LOCATIONS 
C  INTO  AN  IMPULSE  RESPONSE  BY  MULTIPLYING 
C  OUT  THE  NUMERATOR  POLYNOMIAL.  THE  RESULTS 
C  ARE  PRINTED  OUT  AS  POLYNOMIAL  COEFFICIENTS. 

COMPLEX  ZEROI64) 

DIMENSION  ZPOLY (64,3),ZSUM(66),ZTEMP(66) 

K  =  I 

ZSUMI 1 )  =  1. 

DO  90  I  =  2.NZERO+  1 
90  ZSUM(I)  =  0. 

DO  100  I  =  1  ,NZERO 
P=  REAL(ZEROd)) 

Q  =  AIMAG(ZER()(D) 

I!  (Q  I  T  0.)  GO  TO  100 


20  1)0  38  I  —  1  ,N 
w=2*pr<i-n/N 

/-COS(W)  +  J’SIN(W) 

IINUM  =  ( 1,0  ) 

II (N/.ERO.EQ.O)  Ci()  TO  34 
DO  33  K=  l,NZERO 

33  HNUM  =  HNUM*(Z-ZERO(K)) 
HNUMMAG=SQRT(REAL(HNUM)**2  + AIMAG(HNUM)**2) 
IF  (HNUMM  AG.GT.HM  AX)  HMAX=HNUMMAG 

34  HDEN  =  (I.,0.) 

IF(NPOLE.EQ.O)  GO  TO  38 
DO  36  K  =  l.NPOLE 

36  I  (DEN  =  1IDEN*(Z-P0LE(K)) 

38  H(I)  =  HNUM/HDEN 
DO  39  1  =  1, N 

39  11(1)  =  H(I)/HMAX 
RETURN 

END 

■■>t^^*t***************************  ******  ******* 


SUBROUTINE  COEF(PCOEF,NPOLE,ZCOEF\NZERO,N,H) 
C  THIS  ROUTINE  CALCULATES  THE  FREQUENCY 
C  RESPONSE,  H(W),  FROM  THE  TRANSFER  FUNCTION 
C  NUMERATOR  AND  DENOMINATOR  COEFFICIENTS. 
COMPLEX  J,Z 

COMPLEX  HNUM, HDEN,  11(1024) 

DIMENSION  PCOEF (64) ,ZCOEF(64) 

Pf  =  4*ATAN(l.) 

J  =  (0,1) 

C=l. 

IF  (NPOLE.EQ.O)  GO  TO  21 
DO  20  K  =  2,NPOLE  +  I 

20  C  =  C'-ABS(PCOEF(K)) 

21  DO  28  1  =  1, N 
W  =  2*PI*(1-1  )/N 
Z=COS<W)+J*SIN(W) 

IF  (NZERO.GT.O)  GO  TO  22 
IINUM  =  ( I  ,0.) 

GO  TO  24 

22  HNUM  =  (ft.,0.) 

DO  23  K  =  I  ,NZERO+  I 

23  HNUM  =  HNUM  +  ZCOEF(K)*Z**(K-l ) 

24  IF  (NPOLE.GT.O)  GO  TO  25 
IIDHN  =  ( 1,0.) 

CiO  TO  28 

25  IIDEN  =  (0  ,0.) 

DO  26  K  =  LNPOLF+  I 

26  IIDI;.N  =  IIDEN  +  PC()EF(K)*Z**(K-I) 

28  H(I)  =  C*IINUM/HDEN 

RETURN 

END 


C 


( 

SUBROUTINE  M AG(X,Y,Z,PHI,N) 
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ZPOLYIK.1 )  =P**2  +  Q**2 
ZP0LY(K,2)  =  -2.*P 
ZPOLY  (K,3)  =  1 , 

K  =  K  +  1 
100  CONTINUE 
DO  130  1  =  1.  K-l 

ZTEMP(  1  >  =ZSUM(1  )*ZPOLY(l,l ) 

ZTEMP(2)  =ZSUM(2)*ZPOLY  (1,1 )  +  ZSUM(  1  )*ZPOLY(l,2) 
DO  1 10  J  =  3.2*1  -H  1 

110  ZTEMP(J)  =ZSUM(J)*ZPOLY(l,  I )  +  ZSUM(J-1  )*ZPOLY(I,2) 
+  ZSUM(J-2)*ZPOLY(I,3) 

DO  120  L  =  l. 2*1  +  1 
120  ZSUM(L)  =  ZTEMP(L) 

130  CONTINUE 
print  *,  '  ' 

print  \  'POLYNOMIAL  COEFIC1ENTS:’ 

DO  30  1  =  1.NZERO+  1 
30  WRITE  (6,40)  I-l.ZSUMU) 

40  FORMAT  (1X,’Z**’,I2,’  COEF  =  ’,F9.5) 
print  *,  '  ' 

RETURN 

END 

C 

^****** ************************************** 

C 

SUBROUTINE  FREQPZ(ZFREQ,NFREQ,ZERO,NZERO) 

C  THIS  ROUTINE  CONVERTS  ZERO  LOCATIONS 
C  STATED  AS  FRACTIONS  OF  THE  CLOCK  FREQUENCY 
C  TO  Z-PLANE  ZERO  LOCATIONS 
C  NOTE:  THEY  MUST  BE  ON  THE  UNIT  CIRCLE 
COMPLEX  ZERO(64) 

DIMENSION  ZFREQ(64) 

NZERO  =  0 
PI  =4*  ATAN(  1 .) 

DO  10  I  =  l.NFREQ 
W  =  2*PI*ZFREQ(I) 

NZERO  =  NZERO+  1 

ZERO(NZERO)  =  CMPLX(COS(W),SIN(W)) 

IF  ( ABS(S1N(W)).LT.1  E-6)  GO  TO  10 
NZERO  =  NZERO+  1 

ZERO(NZERO)  =CMPLX(COS(W),-SIN(W)) 

10  CONTINUE 
RETURN 
END 
C 

£ ************************************** 

SUBROUTINE  PZ(POLE,NPOLE,ZERO,NZERO,N,H) 

C  THIS  ROUTINE  CALCULATES  THE  FREQUENCY 
C  RESPONSE,  H(W),  FROM  THE  POLE  AND  ZERO 
C  LOCATIONS. 

COMPLEX  J.Z 

COMPLEX  HNUM,HDEN,H(I024) 

COMPLEX  POLE(64).ZERO<64) 

PI  “4*A  I  AN(  I .) 

J  =  (0.,l.) 

UMAX  =0. 


I  I  k 
k  k/* 

(»<)  1 1  Mi 
7  J  -  J  4  k 
PI  =4*  ATAN  ( 1 .) 

DO  30  L=I,M 
LE  =  2**L 
LEI  =  LE/2 
U  =  (l.  0,0.0) 

W  =  CMPLX(C0S(P1/FL0AT(LE1)),S1N(P1/FL0AT(LE1))) 
DO  20  J  =  1  ,LE1 
DO  10  1  =  J,N,LE 
IP  =  I  +  LE1 
T  =  X(IP)*U 
X(IP)  =  X(I)-T 
10  X(I)  =  X(I)+T 
20  U  =  U*W 
30  CONTINUE 
RETURN 
END 

******************************************* 


SUBROUTINE  PZMINM  AX(POLE,NPOLE,ZERO,NZERO,PCOEF,N,Y) 
C  THIS  ROUTINE  PRINTS  THE  PEAKS  AND  VALLEYS  OF 
C  THE  FREQUENCY  RESPONSE  FROM  THE  POLE  AND  ZERO 
C  LOCATIONS 

DIMENSION  Y ( 1 024 ) 

DIMENSION  PCOEF(64) 

COMPLEX  J,Z 

COMPLEX  XNUM,XDEN,X 

COMPLEX  POLE(64),ZERO(64) 

IRANCiE  =  4 
PI=4*ATAN(1.) 

print  \  ’  FREQ  ANGLE  MAG  MAGDB' 

C 

C  CALCULATE  NORMALIZING  COEFFICIENT  ’C' 

C 

C  =  1 . 

IF  (NPOLE.E Q.0)  GO  TO  8 
DO  6  K  =  2,NPOLE+  1 
6  C  =  C-ABS(PCOEF(K)) 

C 

C  IF  MAG  IS  INCREASING  L=1 
C  IF  MAG  IS  DECREASING  L  =  0 
(  INITIALLIZE  L  =  2 
C 

8  L  =  2 

DO  90  1  =  2, N/2 
A  =  Y(l  +  I ) 

B  =  Y(l) 

IF  ( A.GT.B)  GO  TO  10 
IF  (L.EQ.OI  GO  TO  90 
L=0 

print  *,  'MAX' 

GO  TO  20 
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C  THIS  ROUTINE  CALCULATES  THE  ABSOLUTE 
C  MAGNITUDE,  MAGNITUDE  IN  DB,  AND  PHASE  OF 
C  A  COMPLEX  ARRAY  (EG.  FREQUENCY  RESPONSE  H(W)). 
COMPLEX  XG024) 

DIMENSION  Y(1024),/(  1024), PHK 1024) 

PI  =  4*  AT  AN  ( 1.) 

DO  50  1  =  1, N 
XRE  =  REAL(X(D) 

XIM  =  AlMAG(Xd) ) 

Y(I)  =  SQRT(XRE**2  +  XIM**2) 

PHI(I)  =  (180./P1)*ATAN(X1M/XRE) 

IFUXRE.LT.O.).  ANDTXIM.LT.O.))  PHI(I)  =  PHKD-180. 
IF((XRE.LT.O.).AND.(XIM.GT.O.))  PHK1)  =  PHI(1)  +  180. 
1F(ABS(Y ( I ) ) . LT.  1 E-6 )  Y(l)  =  !E-6 
50  Z(I)  =  20*ALOG10(Y(D) 

RETURN 
END 

****************************************** 


SUBROUTINE  SAMPLE(Y,Z,PIII,N) 

THIS  ROUTINE  MODIFIES  THE  FREQUENCY 
RESPONSE  BY  SIN  X  /  X  FOR  SAMPLING  AT  THE 
CLOCK  FREQUENCY 
DIMENSION  Y ( 1 024) ,Z(  1 024) , PHI ( 1 024) 

Pl  =  4*ATAN(l.) 

DO  51  I  =  2,N/2  +  1 
XJ  =  (I-I)*PI/N 

Y(I)  =  Y(l)*(  1-0.8*  (1 -SIN  (XJ)/XJ))) 

51  Z(I)  =  20‘ALOG10(Y(I)) 

RETURN 

END 

******************************************** 


SUBROUTINE  FFT(X,N) 

C  THIS  ROUTINE  PERFORMS  AN  N-POINT  EFT  ON 
C  THE  COMPLEX  ARRAY  X.  AND  RETURNS  THE 
C  RESULT  IN  X. 

COMPLEX  X ( I024),U,W,T 
M  =0 
NN  =  N 

DO  WHILE  (NN.GT.I) 

NN  =  NN/2 
M  =  M  +  1 
END  DO 
NV2  =  N/2 
NM1  =N-I 
J  =  1 

DO  7  1  =  1, NM1 
IF(I.GE.J)  GO  TO  5 
T  =  X  ( J ) 

X(JI  =  X(I) 

X  (1 )  =T 

5  K  =  NV2 

6  II  (K.GE.J)  CiO  TO  7 
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m  II  (I  I  q  1 1  go  in  >»(i 

I  -  l 

prim  *,  'MIN' 

C 

C  PRINT  OUT  MAX  OR  MIN  VALUES 
C  AND  VALUES  AROUND  MAX  OR  MIN 
C 

20  WW  =  2*PI*(I-1  )/N 
DO  80  IW  = -I  RANGE, I  RANGE 
W  =  WW  +  (2*PI*1W)/ (N*l  RANGE) 

Z  =  COS(W)  +  (0.,1.)*SIN(W> 

XNUM  =  ( I  .,0.) 

II  (NZERO.EQ.O)  GO  TO  40 
DO  30  K=  l,NZERO 

30  XNUM  =  XNUM*(Z-ZERO(K)) 

40  XDEN  =  (1.,0.) 

II  (NPOLE.EQ.O)  GO  TO  70 
DO  60  K=  1,N POLE 
60  XDEN  =  XDEN*(Z-POLE(K) ) 

70  X=C*XNUM/XDEN 
XRE  =  REAL(X) 

XIM  =  AIM  AG(X) 

XM  AG  =  SQRT(XRE**2  +  XIM**2) 

XMAGDB  =  20*ALOCil0(XMAG) 

80  WRITE  (6,85)  W/(2*PI),W*360/(2*PI),XMAG,XMAGDB 
85  FORM AT(E13,4,F9.2,E12.3,F10.3) 

90  CONTINUE 
RETURN 
END 
C 

ft******************************************* 

( 

SUBROUTINE  PRNTIMP(T,NT) 

C  THIS  ROUTINE  PRINTS  THE  IMPULSE  RESPONSE. 
DIMENSION  T ( 1 024) 
print  *,  ' 

print  *.  'IMPULSE  RESPONSE:’ 
print  \ 

WRITE  (6,15)  <I,T(I),I  =  I,NT) 

15  FORMAT(l5,l  l2.5) 
print  *, 

RETURN 

END 

( 

£•**.  ****.**«.«**«****««*«*«»***•***»**»«*•** 

C 

SUBROUTINE  PRNTI  RI  Q(Y,Z,PHI,N) 

(  THIS  ROUTINE  PRINTS  THE  FREQUENCY  RESPONSE. 
DIMENSION  Y (1024), Z(  1024), PIIK 1024) 
print  *,  '  ' 

print  *,  FREQUENCY  RESPONSE:' 
print  *,  '  ’ 

WRITE  (6,24) 

24  FORMAT  (3X,'#',I0X,’M AG’,8X,’M AGDB’,8X,’PH ASE’) 
WRITE  (6,25)  (l,Y(l),Z(l),PHI(l)J  =  UN/2  +  l) 

25  FORMAT  (!5,3X,FI 2.5,3X,F10.3,3X,FI0.3) 
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print  *, 

RETURN 

END 

C 

^#***,.*#*****, ****************************** 

C 

SUBROUTINE  PRNTPZIPOLE.NPOLE.ZERO.NZ.ERO) 

C  THIS  ROUTINE  PRINTS  THE  POLE  AND  ZERO 
C  LOCATIONS  IN  X-Y  AND  R-THETA  COORDINANTS. 

COMPLEX  POLE(64),ZERO(64) 

PI  =  4*  AT  AN  (1.) 

IE  (NPOLE.EQ.O)  GO  TO  30 
print  *. 

print  *,  POLE  LOCATIONS’ 

print  V  #  REAL  IMAG  RADIUS  ANGLE  FREQ  (EC)' 
DO  25  I  =  l.NPOLE 

RADIUS  =  SQRT(REAL(POLE(l))**2  +  AIM  AG(POLE(l))**2) 

P  =  REAL(POLEd)) 

Q  =  AIM  AG  (POLE  (I )) 

IE  (ABS(P).GT.lE-8)  GO  TO  15 
ANGLE  =  90. 

IE  (Q.LT.O)  ANGLE  =  -90. 

GO  TO  20 

15  ANGLE  =  ( 180./PD*  ATAN(  AIMAG(POLEd))/REAL(POLE(l))) 

IF  (( AIM  AG(POLE(1)).GT.O.).  AND.(REAL(POLE(1)).LT.O.)) 
&ANGLE  =  ANGLE  + 180. 

IE  ( ( AIM  AG(POLE(I)).LT.O).AND.(REAL(POLE(1)).LT.O.)) 
&ANGLE  =  ANGLE- 180. 

IF  ( ANGLE. LT. -1 E-6)  GO  TO  25 
20  WRITE(6,10)  l.POLE(l), RADIUS, ANGLE, ANGLE/360. 

10  FORM  AT  (1 5,5  FI  0.4) 

25  CONTINUE 

30  IF  (NZERO.EQ.O)  RETURN 
print  *,  ' 

print  *.  ZERO  LOCATIONS’ 

print  #  REAL  IMAG  RADIUS  ANGLE  FREQ  (EC1)' 
DO  45  I  =  UNZERO 

R  ADIUS  =  SQRT(REAL(ZERO(l))**2  +  AIM  AG(ZERO(l))**2) 

P  =  REAL(ZERO(D) 

Q  =  AIMAG(ZEROd)) 

IF  <ABS(P).GT.lE-8)  GO  TO  35 
ANGLE  =  90. 

IF(Q.LT.O)  ANGLE  =  -90. 

GO  TO  40 

35  ANGLE  =  (180./Pi)*ATAN(AIMAG(ZERO(I))/REAL(ZEROd))) 

IE(  ( AlMAG(ZEROd)).GT.O).  AND.(REAL(ZERO(I)).LT.O.)) 

&  ANGLE  =  ANGLE  +180 

II  ( ( AIM ACi(ZERO(l)).LT.O.).AND.(REAL(ZERO(I)).LT.O.)) 

&  ANGLE  =  ANGLE- 1 80. 

II  (ANGLE  LT.-IF.-6)  GO  TO  45 
40WRHI  ((>.10)  l./I  ROd), RADIUS  ANtill  .ANGI.1 73(>0 
I  >  <  <  IN  I  |NU| 

I'Miii 

Ul  I  URN 
I  Nl) 


A  2- 7 


o  n 


Appendix  3 

BANDPASS  FILTER  DESIGN  PROGRAM 


c  CDF4  5/26/83 
c 

c  This  program  plots  the  frequency  response  for  cdf4 

c 

COMMON  KQUEST,KPRNT,KLAB,SCALE 
COMMON  Y ( 1024), Z(  1024),PHI0024),N,H 
COMMON  /ZZZ/PCOEF(64),NPOLE,ZFREQ(64),NFREQ 
COMPLEX  G(1024),H(1024) 

DIMENSION  XPOINT(64),YPOINT(64) 

PI=4*ATAN(1.) 

M=  10 
N  =  2**M 
SCALE  =  -1.0 
KPRNT  =  0 
KQLJEST  =  1 
C 

+**+*******+*****++**++**+***+++***++*+++4+**+ 

print  *,  '  ’ 
print  *,  'TITLES’ 

DO  8  1  =  1,  N 

8  G(l)  =  (1,0.) 

KLAB  =  100 

CALL  PLTPZ(POLE,0,ZERO,0,KLAB,SCALE) 

CALL  PLTMG  (Z,0,KLAB,SCALE) 

C 

£•**4 t****************************************** 

C 

10  print  *,  ’  ’ 

print  *,  '91  DEG  POLES' 

NPOLE= 1 7 
DO  12  1  =  1.NPOLE+  1 
12  PCOEF(I)  =0. 

PCOEFO )  =  1 . 

PCOEF(15)  =  .6 
PCOEFO  8)  =3 
KLAB  =  30 
CALL  CALCPOLE 
DO  14  1  =  1, N 
14  G  ( I )  =G(I)*H(I) 

C 

£*+***  ******** ************************************ 

C 

20  print 

print  *,  '89  DEG  POLES’ 

NPOLE=  19 
DO  22  I  =  1  .NPOLE+  1 
22  PCOEFtl)  =0. 

P<  Oil  (I  (  I 

p<  ori  <  i  s»  t, 

PCOl  I  06>=  1 
PCOl.  F 1 20  >  =  2 


KLAB  =  40 
CALL  CALC"  OLE 
DO  24  1=  l,N 
24  CJ  (I)  =  (ti  (l)-ll(D) 


30  print  *,  '  ’ 

print  *,  '90  DEG  POLES' 
NPOLE  =  2 
DO  32  1=1, NPOLE 
32  PCOEF(l)  =0. 

PCOEF(  1 )  =  1 . 

PC’OEF(3)  =  875 
KLAB  =  20 
CALL  CALCPOLE 
DO  34  1=1, N 
34  G(l)=G(l)*II(l)’H(l) 


40  print  *,  '  ' 

print  *,  ZEROES’ 

NP()LE  =  0 
N  FREQ  =  6 

DATA  (ZFREQ(I),  1=1,6)  /.032S,.  108,. 1 80,  320,  392,  4675/ 
KL  AB  =  20 
CALL  CALCZERO 
DO  46  I  =  I ,  N 
46  GO)  =G(I)*H(1) 


50  print  *, 

print  *,  'TOTAL  RESPONSE’ 

KLAB=  10 

IF  (KQUEST.EQ.O)  KLAB  =  KL AB+  1000 
CALL  M AG(G,Y,Z,PHI,N) 

CALL  PLTMG(Z,N,KLAB,SC'ALE) 

60  STOP 


SUBROUTINE  PLTMG(Z,N,KLAB,SCALE) 

DIMENSION  Z(  1024) 
print  *, 

print  \  FREQUENCY  RESPONSE  PLOT’ 

(  ALL  PLT(Z,N,KLAB,SCALF.,0.,.5,.  1  ,-80.,0.,20.,0.,0.,0,0) 
RETURN 


SUBROUTINE  PLTPH(PHI,N,KLAB,SCALE) 
DIMENSION  PIIK1024) 


print  *,  '  ' 

print  \  ’PHASE  RESPONSE  PLOT’ 

CALL  PLT(PHI,N,-KL  AB,SCALE,.22,.28,.01 ,- 1 80  , 180., 90.. ().,()., 0,0) 

RETURN 

END 

C 

c 

SUBROUTINE  CALCPOLE 
COMMON  KQUEST,KPRNT,KLAB,SCALE 
COMMON  Y ( 1024),Z(1024),PHI(  1024),N,H 
COMMON  /ZZZ/PCOEF(64),NPOLE,ZFREQ(64),NFREQ 
COMPLEX  H(1024),POLE(64),ZERO(64) 

NZERO  =  0 

IF  (KQUEST.EQ.O)  KLAB  =  KLAB+ 1000 
CALL  ZPOLR(PCOEF,NPOLE,POLE,lER) 

CALL  COEF(PCOEF,NPOLE,PCOEF,0,N,H) 

IF  (KPRNT.EQ.I)  CALL  PRNTPZ(POLE,NPOLE.ZF.RO,NZERO) 

CALL  PLTPZ(  POLE, NPOLE,ZERO,NZERO,KLAB,SC  ALE) 

CALL  MAG(H,Y,Z,PHI,N) 

CALL  PLTMC.(Z,N,KLAB,SCALE) 

IF  (KPRNT.EQ.2)  CALL  PZMINM AX(POLE,NPOLE,ZERO,NZHRO,PCOEF,N,Y 
RETURN- 
END 
C 

** * ******** * * * * * * * * ** * ******** * **** * *** * 

C 

SUBROUTINE  CALCZERO 
COMMON  KQUEST,KPRNT,KLAB,SCALE 
COMMON  Y ( 1024), Z(  1024), PHI (1024), N,H 
COMMON  /ZZZ/PCOEF(64),NPOLE,ZFREQ(64),NFREO 
COMPLEX  HU024),POLE(64),ZERO(64) 

NPOLE=0 

IF  (KQUEST.EQ.O)  KLAB-KLAB+ 1000 
CALL  FREQPZ(ZFREQ,NFREQ,ZERO,NZERO) 

CALL  PZ(POL.E,NPOLE,ZERO,NZERO,N,H) 

IF  (KPRNT.EQ.I )  CALL  PRNTPZ(POLE,NPOLE,ZERO,NZERO) 

CALL  PLTPZ(POLE,NPOLE,ZERO,NZERO,KLAB,SCALE) 

CALL  MAG(H,Y,Z.PHI,N) 

CALL  PLTMG(Z,N,KLAB,SCALE) 

II  (KPRNT.EQ.2)  CALL  PZMINMAX(POLE,NPOLF.,ZERO,NZERO,N,Y) 

IF  (KPRNT.EQ.I)  CALL  PZ!MP(ZERO,NZERO) 

RETURN 

END 


Appendix  4 
WHITE  PAPER 

F\1  DEMODULATION  FOR  FREQUENCY  HOP  COMMUNICATIONS  SYSTEMS 

1.  INTRODUCTION 

for  (he  past  several  years,  (he  General  Electric  Corporate  Research  and  Development  Center  has 
been  developing  Charge  Domain  Filter  (CDF)  technology  for  application  to  Military  Systems  needs 
(Work  performed  under  contract  I- 1  %28-80-C-02 13.)  A  CDF'  is  a  type  of  Charge  Coupled  Device 
(CCD)  for  signal  processing  in  which  samples  of  the  incoming  signal  are  represented  by  charge  packets 
which  are  split  up,  routed  along  various  paths,  and  recombined  to  form  new  charge  packets  representing 
samples  of  the  output  signal.  All  necessary  mathematical  operations  including  multiplication,  addition 
and  delay  can  be  performed,  resulting  in  a  completely  flexible  signal  processing  capability.  Since  all 
operations  are  performed  directly  on  the  charge  packets  themselves,  non-linearities  that  commonly  arise 
in  charge- to- voltage  and  voltage-to-charge  translations  do  not  occur.  These  devices  are  capable  of  im¬ 
plementing  extremely  sophisticated  signal  processing  functions  with  high  accuracy,  low  power,  and  very 
high  speed.  After  working  out  many  details  of  device  structure  design  and  filler  system  architecture, 
this  effort  has  recently  culminated  in  the  design  of  a  bandpass  filter  with  some  unusual  characteristics. 
These  are  illustrated  in  Figs.  I  and  2,  which  are  amplitude  plots  of  the  frequency  response  of  the  filler 
The  width  of  the  passband  is  approximately  1%  of  the  clock  frequency,  and  the  pass  band  has  both  a 
linear  phase  response  and  a  flat  amplitude  response.  The  skirt  selectivity  corresponds  to  8  poles,  and 
the  stop  band  drops  off  to  more  than  -100  dB  compared  to  the  passband  response.  Furthermore,  since 
the  insertion  loss  is  only  about  12  dB,  the  dynamic  range  is  expected  to  be  of  the  order  of  70  dB.  Since 
the  design  is  compatible  with  buried  channel  CCD  structures,  it  is  expected  that  this  filter  characteristic- 
can  be  obtained  at  any  frequency  from  the  audio  range  up  to  approximately  100  Mil/.  Since  the  center 
frequency  of  this  filter  depends  on  the  clock  frequency,  it  is  possible  to  move  its  frequency  around  at 
electronic  speeds.  Thus,  a  filter  similar  to  the  one  shown  could  serve  as  a  basis  for  a  highly  secure 
communications  system  based  on  frequency  hopping. 

2.  POSSIBLE  SYSTEM  ARCHITECTURE 

A  filler  such  as  the  one  shown  can  be  used  in  many  different  applications,  and  the  communications 
system  proposed  here  is  only  one  suggestion.  It  was  chosen  as  an  example  because  it  exploits  a  number 
of  the  unique  characteristics  of  the  CDF  filter  cited  above. 

A  possible  system  architecture  for  a  frequency  hopping  communications  system  is  shown  in  Fig  . la 
and  lb.  Fig. la  is  a  schematic  system  block  diagram  including  a  tuner  and  first  i  f.  amplifier  of  conven¬ 
tional  design  which  delivers  a  relatively  broad  frequency  spectrum  and  a  monolithic  “back  end”  whose 
function  is  to  select  a  frequency  hopped  signal  from  within  the  passband  of  the  front  end  and  deliver  a 
fully  decoded  digital  bit  stream  at  its  output.  Since  the  monolithic  back  end  is  a  sampled  data  system 
(and  therefore  has  alias  passbands),  we  assume  that  the  passband  of  the  front  end  is  limited  to  less  than 
one-fourth  of  the  sampling  frequency  of  the  monolithic  processor.  This  is  indicated  in  Fig. 3b  by  the 
dashed  line  frequency  response.  In  the  monolithic  processor  itself,  a  first  filter  is  used  in  one  of  its 
alias  response  bands  to  select  a  relatively  narrow  portion  of  the  spectrum  and  to  down  convert  it  to  the 
fundamental  passband  of  the  filter  at  one-fourth  of  its  clock  frequency.  This  band  of  frequencies  is 
then  further  split  by  a  second  filter  of  the  same  type  that  is  operated  at  a  lower  clock  frequency.  Since 
this  second  filler  is  also  operated  in  an  alias  passband,  a  further  down  conversion  is  effected,  and  the 
bandwidth  is  further  narrowed  in  the  process.  Finally,  the  charge  packets  arc  put  through  a  charge 
domain  limiter  so  as  to  remove  any  amplitude  modulation  that  may  be  present,  and  is  passed  to  a 
charge  domain  FM  demodulator  which  will  be  discussed  further  below. 

In  the  above  system  diagram,  it  is  assumed  that  all  of  (he  operations  arc  performed  in  the  charge 
domain  so  that  no  signal  recovery  is  required  at  any  internal  point.  The  entire  system  is  assumed  to  be 
implemented  by  splitting  and  recombining  of  charge  packets  and  by  skimming  operations  This  will 


per  mu  a  single  chip  implementation  which  can  operate  at  high  speed  with  good  immunity  to  power  sup¬ 
ple  variations,  etc.  In  order  to  accomplish  this  objective,  several  of  the  component  subsystems,  such  as 
the  Charge  Domain  amplitude  limiter  and  the  Charge  Domain  FM  demodulator  would  need  lurther 
development  effort. 

3.  PROPOSED  RESEARCH  PROGRAM 

Ideas  have  already  been  generated  for  accomplishing  the  objective  outlined  above,  and  these  will 
now  be  described. 

a  Amplitude  limiter. 

A  method  for  removing  amplitude  modulation  from  a  frequency  modulated  signal  is  shown  in 
Fig. 4a  and  4b  Fig.  4a  is  a  surface  potential  diagram  for  a  new  device  structure  which  accomplishes 
amplitude  limiting  of  an  a.c.  signal,  and  Fig. 4b  is  a  plan  view  schematic  layout.  The  operation  is 
based  or,  the  principle  of  skimming. (Ref.l  Charge  packets  representing  samples  of  a  sinusoidally 
varying  input  signal  arrive  from  the  reservoir  at  the  left  of  Fig.4.  and,  assuming  their  amplitude  is 
sufficient  to  surmount  the  skimming  barrier,  they  pass  over  the  small  output  reservoir  on  their  way 
to  the  sink  In  t he  process,  the  output  reservoir  is  filled  completely,  irrespective  of  the  amplitude 
of  the  incoming  packet.  If,  on  the  other  hand,  the  incoming  packet  was  insufficient  to  surmount 
the  skimming  barrier,  no  charge  at  all  will  enter  the  output  reservoir.  Thus,  the  output  reservoir 
will  either  be  filled  completely  or  totally  empty  (except,  of  course  for  a  very  narrow  range  of  ampli¬ 
tudes  very  close  to  the  level  of  the  skimming  barrier).  If  the  level  of  the  skimming  barrier  is  adap¬ 
tively  adjusted  so  that  it  lies  right  at  the  center  of  the  range  representing  the  positive  and  negative 
peaks  of  the  input  signal,  the  result  will  be  to  turn  signals  of  essentialy  any  amplitude  into  square 
waves  of  uniform  amplitude. 

b.  I'M  Demodulator 

The  demodulator  proposed  for  this  system  is  essentially  two  frequency  sensitive  fillers,  one 
tuned  to  the  frequency  excursion  below  the  carrier,  and  the  other  tuned  to  frequency  excursion 
above  the  carrier.  In  order  to  keep  the  transient  response  as  short  as  possible,  these  filters  can  be 
designed  as  transversal  filters  with  impulse  responses  that  are  of  finite  duration.  In  practice,  fillers 
having  impulse  responses  as  short  as  9  samples  are  feasible.  Thus,  digital  bit  rales  can  be  as  high 
as  one  tenth  of  the  clock  frequency  are  possible.  These  filters  can  be  shaped  to  produce  a  linear 
amplitude  versus  frequency  response  in  the  area  of  interest  as  shown  in  Fig. 5.  Or,  for  a  digital  f.m. 
demodulation,  where  only  a  single  frequency  shift  is  present,  these  filters  can  be  more  narrowly 
tuned  and  centered  on  the  two  shifted  frequencies. 

4.  SUMMARY 

Some  excellent  and  unique  capabilities  possessed  by  Charge  Domain  Filters  have  already  been 
demonstrated,  and  it  is  proposed  to  exploit  these  in  a  frequency  shift  (i.e.  FM)  digital  communications 
system.  This  system  would  be  tunable  in  frequency  at  electronic  speeds,  permitting  a  high  speed  fre¬ 
quency  hop  strategy  for  security  and  jam  resistance.  The  entire  "back  end"  including  second  and  third 
i  f  ,  limiter  and  FM  demodulator  can  be  incorporated  onto  a  single  silicon  chip  if  desired. 

The  basic  filter  technology  required  — including  both  device  structures  and  filter  system  architecture 
have  already  been  developed,  and  specific  ideas  for  implementing  the  limiter  and  demodulator  fuctions 
have  been  presented. 


Appendix  5 

PROGRAMMABLE  C  HARGE  DOMAIN  FILTERS  FOR  SOLID  STATE  RADARS 

1.  INTRODUCTION 

During  ihe  past  two  years,  the  General  Electric  Corporate  Research  and  Development  Center  has 
been  developing  Charge  Domain  Device  (CDD)  technology  lor  potential  use  in  military  applications 
These  early  developments,  while  not  addressing  particular  applications,  demonstrated  that  greatly  in¬ 
creased  linearity  and  significantly  higher  clock  speed  (as  compared  to  conventional  techniques)  could  he 
achieved  by  this  new  technology.  The  thrust  of  this  work  was  to  improve  device  operating  speed  ami  ,!t- 
euracy  and  to  investigate  the  underlying  device  physics  which  determines  the  ultimate  performance 
capabilities  of  this  technology.  Device  simulation  tools  have  been  developed  to  study  these  new  device 
structures.  With  the  knowledge  gained  thus  far,  it  is  now  appropriate  to  identify  some  areas  of  applica¬ 
tion  where  CDD’s  can  best  be  used. 

A  series  of  presentations  on  the  subject  of  CDD  capability  was  recently  made  to  groups  of  engineers 
within  GE  who  are  involved  with  military  applications.  These  groups  included  representatives  from 
G.E.'s  Aerospace  Electronic  Systems  Dept.,  Military  Equipment  Systems  Operation,  Ordnance  Systems 
Dept.,  and  other  parties  knowledgeable  in  military  applications  for  sophisticated  electronic  signal  pro¬ 
cessing.  One  frequently  expressed  desire  was  for  flexibility,  and  in  particular,  for  an  ability  to  digitally 
program  the  function  performed  by  a  signal  processing  sub-system.  Our  efforts  have  correspondingly 
included  the  investigation  of  programmable  CDD's,  and  recent  developments  have  indicated  the  feasi¬ 
bility  of  providing  this  important  function  in  CDD  signal  processing  devices.  The  discussions  also 
focussed  on  several  specific  application  areas  which  are  within  the  reach  of  present  CDD  technology, 
but  which  are  extremely  difficult  to  accomplish  by  any  other  means.  A  building  block  approach  has 
been  conceived  which  would  utilize  a  chip-cascadable,  programmable  CDD  analog-analog  correlator  in 
order  to  meet  the  needs  of  many  of  these  programs. 

2.  TECHNICAL  APPROACH 

The  technical  approach  for  designing  a  generic  chip  which  could  meet  these  needs  is  based  on  a  new 
signal  processing  chip  architecture  which  could  be  called  a  Moving  Reference  Programmable  CDD.  In 
this  architecture,  analog  samples  are  scanned  into  a  series  of  individual  programmable  charge  splitters 
where  they  remain  until  they  are  replaced.  Binary  bits  representing  the  K  and  l-K  coefficients  for  a  tap 
weight  propagate  along  digital  shift  registers  so  that  the  coefficients  move  from  splitter  to  splitter  in  se¬ 
quence.  During  each  bit  time,  the  charge  packet  in  each  of  the  splitters  is  split  into  two  equal  parts, 
and  the  portion  corresponding  to  the  binary  weight  of  that  bit  is  delivered  to  the  K  or  (l-K)  accumula¬ 
tor  as  appropriate.  The  remaining  half  is  then  ready  for  further  binary  splitting.  When  all  of  the  binary 
bits  of  a  particular  tap  coefficient  have  been  processed,  the  output  of  the  entire  scries  is  read  out  by 
means  of  a  pair  of  overlying  electrodes  that  sense  the  charges  in  the  K  and  l-K  accumulators.  Alter 
readout,  the  charges  in  the  K  and  l-K  accumulators  are  returned  to  the  same  charge  splitting  reservoir 
where  they  started.  That  is,  they  are  not  passed  to  the  next  stage  as  they  would  be  in  a  conventional 
CCD  filter.  At  this  point,  the  shift  register  containing  the  reference  coefficients  is  clocked  so  that  the 
binary  bits  move  to  the  next  splitter  in  the  series.  The  advantages  of  this  architecture  include  the  fact 
that  there  is  no  cumulative  charge  transfer  loss,  the  digital  bits  can  be  passed  from  chip  to  chip  without 
degradation,  and  the  output  can  span  several  chips,  if  desired.  The  lack  of  charge  transfer  loss  is  due  to 
the  fact  that  the  analog  charge  packets  do  not  move  from  stage  to  stage.  Thus,  any  charge  left  behind  as 
a  result  of  charge  transfer  inefficiency  is  not  mixed  with  the  next  sample,  but  is  instead  simply  recom¬ 
bined  with  the  transferred  pari  of  the  same  sample.  The  cascading  of  chips,  which  is  required  for  the 
long  BT  products  needed  in  many  applications,  is  accomplished  by  simply  connecting  the  shift  register 
outputs  from  one  chip  to  the  inputs  of  (he  next  anti  connecting  Ihe  corresponding  overlying  electrodes 
together. 


Assuming  that  H  to  Id  bits  <>l  .uuii.kv  ale  needed  on  ilic  relercncc  coefficients  lor  adequate 
side  lube  suppression,  ihe  splm my.  clock  will  have  to  run  at  8  to  10  limes  the  primary  sampling  rale. 
This  is  not  a  problem  however,  since  the  splitting  clock  on-chip  driver  would  only  be  driving  one  low- 
capacitance  electrode.  Additional  speed  may  be  obtained  by  implementing  the  splitting  function  in 
buried  channel  while  the  output  charge  sensing  procedure  (which  runs  at  the  sampling  rate)  could 
remain  surface  channel  for  linearity  reasons.  Thus  sampling  rates  in  the  multi-megahertz  range  are  en¬ 
tirely  feasible. 

.1.  DISCUSSION 

The  development  of  such  a  correlator  chip  would  provide  a  practical  solution  to  a  broad  range  of 
problems  including  a  programmable  sampled  data  filter,  a  pulse  compressor  for  solid  state  radars,  and  a 
secure  communications  device,  incorporating  coded  variable  analog  reference  techniques  as  well  as  fre¬ 
quency  hopping  techniques.  One  particular  application  which  is  impractical  with  conventional  technolo¬ 
gy  is  outlined  below. 

At  the  present  time,  USAF  advanced  warning  airborne  radars  utilize  high  power  pulsed  klystrons 
with  peak  power  in  the  multi-megawatt  range  in  order  to  provide  sensitivity  at  the  desired  distances. 
Not  only  are  these  lubes  bulky  and  relatively  heavy,  but  they  are  also  very  expensive  and  difficult  to 
build  A  desirable  alternative  would  be  an  equivalent  solid  state  radar  power  source,  but  a  direct  re¬ 
placement  is  out  of  the  question  because  of  the  limitations  on  peak  power  of  solid  state  devices.  In¬ 
stead,  the  solid  stale  alternative  would  use  a  pulse  of  longer  duration  and  lower  peak  power,  with  ap¬ 
propriate  pulse  compression  to  obtain  the  equivalent  peak  power.  The  problem  comes  about  because  of 
the  extreme  amount  of  pulse  compression  required,  (namely  several  thousand  to  one),  the  high 
bandwidth,  (namely  several  MHz  ),  and  the  low  sidclobe  level  that  can  be  tolerated.  These  require¬ 
ments  in  combination  add  up  to  several  giga  multiply/add  operations  per  second  with  at  least  12  bits  ac¬ 
curacy.  Unfortunately,  the  size  and  weight  of  a  processor  to  accomplish  this  signal  processing  task  is 
presently  at  least  as  large  as  the  klystron  it  would  replace,  and  the  power  required  would  be  consider¬ 
ably  greater  than  the  transmitter  power.  By  using  programmable  and  cascadable  CDD  correlators,  this 
task  can  be  accomplished  with  perhaps  a  few  dozen  low  power  chips. 

The  alternative  of  addressing  this  application  with  conventional  CCD's  is  not  an  attractive  one.  If 
surlace  channel  CCD's  were  used,  charge  transfer  losses  would  be  prohibitive.  Even  buried  channel 
devices  become  marginal  after  so  many  stages,  but  these  could  be  made  to  work  if  a  dummy  stage  (with 
bias  charge  only)  were  inserted  between  each  signal  sample.  This  means  that  the  clock  rale  would  have 
to  be  double  the  sample  rate,  implying  clock  rales  as  high  as  10  MHz.  in  some  cases.  The  number  of 
stages  is  also  doubled,  of  course,  so  the  power  is  increased  fourfold.  We  would  end  up  with  a  CCD 
with  about  5000  stages  (10,000  transfers)  being  clocked  at  10  MHz.  To  avoid  transfer  loss  problems,  a 
CCD  of  this  length  and  speed  would  have  to  be  a  buried  channel  device,  but  this  would  imply  a  high 
degree  of  output  non-linearity.  In  other  words,  this  application  is  not  well  suited  to  any  previously  pro¬ 
posed  technology -either  digital  or  analog. 

4.  CONCLUSION 

Ihe  development  of  Charge  Domain  Devices  has  reached  the  stage  where  they  may  be  used  in  ap¬ 
plications  to  overcome  problems  that  are  extremely  cumbersome,  if  not  insoluble,  with  conventional 
techniques  One  particular  application  which  exemplifies  this  is  pulse  compression  for  solid  stale  long 
range  radars.  The  high  speed  and  programmability  capabilities  of  C'DD’s  provide  an  attractive  solution 
to  this  problem.  A  follow-on  program  of  approximately  18  months  total  duration  to  develop  this  appli¬ 
cation  would  be  appropriate  at  this  time. 
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